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1.  Introduction 

Several  aspects  of  sensitivity  analysis  in  nonlinear  programming 
have  been  examined  from  a computational  viewpoint  by  Arm&coet  and  Flacco 
(197*0.  That  work  was  based  on  the  theory  developed  by  Flacco  (1973)  and 
used  the  computational  procedures  implemented  by  Armacost  and  My lander 
(1973)  using  the  SUMT-Version  k computer  code  with  the  logarithmic- 
quadratic  loss  penalty  function  modified  to  estimate  the  partial  derivatives 
of  the  solution  point  and  the  objective  function  taken  with  respect  to 
certain  specified  problem  parameters,  Flacco  (1973)  developed  the  necessary 
formulas  to  provide  estimates  of  the  partial  derivatives  of  the  Lagrange 
multipliers  taken  with  respect  to  the  problem  parameters  using  the 
logarithmic-quadratic  loss  penalty  function,  Armacost  and  Flacco  (1975) 
developed  the  formulas  to  obtain  first  and  second  order  sensitivity  infor- 
mation for  the  optimal  value  function  (a  function  of  the  parameters 
defined  by  the  objective  function  evaluated  at  the  solution  point). 
Additionally,  Armacost  and  Flacco  (1976)  have  shown  that  when  the  parameters 
are  the  right-hand  side  components  of  the  constraints,  the  partial  deriva- 
tives of  the  Lagrange  multiplier*  are  the  components  of  the  Hessian  of  the 
optimal  value  function.  The  supporting  tlieory  is  addressed  In  Section  2, 
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Amacont  am  Fiacco  (1^7*0  focus  on  the  computational  experience 
with  four  example  problems  discussing  various  aspects  of  the  sensitivity 
analysis  procedure  and  results.  They  note  that  for  large  problems  with 
a large  number  of  parameters,  a very  large  number  of  partial  derivatives 
will  be  estimated.  This  is  not  only  time  consuming,  but  is  also  burden- 
some to  the  user  who  must  evaluate  all  of  them,  many  of  which  are  zero 
or  very  close  to  zero  in  value.  In  addition,  it  is  unlikely  that  the 
value  of  the  solution  (the  value  of  the  objective  function  evaluated  at 
the  solution  point)  will  be  sensitive  to  more  than  a relatively  few  para- 
meters. Because  of  this,  the  method  developed  by  Armacost  and  Fiacco  (1973) 
to  estimate  the  firrt  order  sensitivity  of  the  optimal  value  function  is 
Incorporated  In  the  computer  program  here  to  provide  an  option  for  pre- 
liminary screening  of  the  parameters  to  determine  which  ones  affect  the 
optimal  value  function.  Using  the  formulas  developed  by  Fiacco  (1973)* 
a second  option  Is  Included  which  permits  the  calculation  of  the  sensitivity 
estimates  for  the  Lagrange  multipliers.  The  computer  code  and  options  used 
to  accomplish  these  calculations  are  discussed  In  Section  3. 

In  Section  **,  the  new  computational  experience  with  Lagrange 
multiplier  and  optimal  value  function  sensitivity  using  these  options  is 
presented  for  the  four  examples  of  Armacost  and  Fiacco  (197**)# 

In  Section  5 • a sensitivity  analysis  Is  conducted  for  a large 
scale,  multi-item  Inventory  model  developed  by  Schrady  and  Choe  (1971) 
for  the  U.  S,  Navy.  Vhlle  the  example  used  Is  the  same  small  one  used  by 
Schrady  and  Choe,  it  nonetheless  exhibits  the  value  of  performing  such  a 
sensitivity  analysis  In  real  world  situations  and  Illustrates  the  care 
that  must  be  taken  In  Interpreting  the  sensitivity  results. 
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2,  Supporting  Theory 

The  problems  considered  here  are  of  the  form  of  Problem  p(C), 
minimize  f(x,€) 

subject  to  g^x,*:)  > 0,  i*l,  • • • ,ro,  r(C) 

h.(x,€)  - 0,  j-1 

J 

When  certain  assumptions  are  satisfied,  Fiacco  (1973)  and  A mac os t 
and  Fiacco  (1975)  have  shown  the  existence  of  the  first  order  sensitivity 
of  a Kuhn-Tucker  triple  and  the  first  and  second  order  sensitivity  of  the 
optimal  value  function.  Additionally,  they  provide  the  means  of  estimating 
this  sensitivity  by  way  of  the  logarithmic-quadratic  loss  penalty  function. 
The  following  four  assumptions  are  sufficient  to  establish  these  results 
and  are  assumed  to  hold  throughout  this  Section, 

Al  — The  functions  defining  Problem  P(€)  are  twice  continuously 
differentiable  in  (x,f)  in  a neighborhood  of  (x*,0). 

A2  —The  second  order  sufficient  conditions  for  a local  minimum 
of  Problem  p(o)  bold  at  x*  with  associated  Lagrange  multi- 
pliers u*  and  w*. 

A3  —The  gradients  for  all  i auch  that  ^(x^O)  * 0, 

and  7xh^(x*,0),  >l,...,p  are  linearly  independent, 

A**  —Strict  complementary  slackness  holds  at  (x*,0)  (i,e.,  u^*  >0 
for  all  1 such  that  ^(x^O)  - 0), 

The  main  results  are  presented  without  proof  and  are  stated  latre 
for  completeness.  The  portion  of  the  theory  used  in  the  computational 
alrorlthn  is  made  specific.  The  Ugrangian  for  iroblem  1(0  is 

B P 

L(x,u,-,€)  • - £ u.c,(x,<)  ♦ £ x,h,(x,*) 

l-l  11  j'l  J ' 

w}*re  u.,  i*i,...  m and  w.#  j»i,...,pare  the  Lagrange  multi  pliers 
- J 
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assoc iated  with  the  Inequality  and  equality  constraints  respectively. 

The  first  result  was  proved  as  Theorem  2,1  by  Fiacco  (1V73)* 
THEOREM  1*  (First  order  sensitivity  of  a Kuhn-Tucker  triple) 

If  assumptions  At,  A3  and  A^  hold  for  Problem  p(f)  at  (x*,0),  then 
(a)  x*  is  a local  isolated  minimizing  point  of  Problem  P(0)  and  the 
associated  I-af'mnrc  multipliers  u*  and  vr*  are  uniquei 
(t)  for  f in  a noirhbortiood  of  0,  there  exists  a unique,  once  con- 

tiriously  differentiable  vector  function  y(()  - (x(f) ,u(C) ,w(f )) 
rutisfyinr  the  second  order  sufficient  conditions  for  a local 
minimum  of  iTotlen  r( ' ) such  that  (x(0)  ,u(0)  ,w(0))  * (x*,u*,w*) 
and  fence,  x(0  is  a locally  unique,  local  minimum  of  Iroblem  J<€) 
with  associated  unique  La*Traur.e  multipliers  u(£)  and  x(C)j  and 
(c)  for  f near  , the  set  of  nindint'  inequalities  is  unchanged, 

strict  complementary  slackness  ijolds  for  u^(C)  for  i such  that 
■^(x(0#0  * 0,  and  the  bindl.-v  constraint  gradients  are  linearly 
independent  at  x(f), 

bet  y(0  » (y(C;,u(-').w(r))  **  a Kuhn-Tucker  triple  where  x(f) 

solver,  i rol-le^  HO,  then  the  optimal  value  function  is  defined  as  f#(C)  ■ 
r[y.( ( \ and  the  ordinal  value  Uu;  ran.*  Ian  Is  l.*(0  • L(x(C)*u(0#w(0»0* 
l.e  secofid  result  was  recently  established  by  Armacost  and  Klacco 
(V>7j;  In  their  Theorem  , stated  here  an  Theorem  2. 

THEOREM  2 t (First  and  second  order  cttfinces  in  the  optimal  value  function) 
If  assumption.  Ai,  A2,  A3  and  A^  l>old  for  Problem  P(0  (x*,0),  then 

for  ( near  0, 


(a)  f(0  - l*(0* 

(b)  7,f*(*  ) * 7rl(x,u,w/)| 

|(x,u,w)Hx(f)iu(r)iw(e)) 


T-335 


m 

" - z u V Wx,€) 

i=l 
P 

+ L w 7 h (xf€)|  ; 

>1  3 3 l(xfulw)-(x(€)fu(f)tw(C)) 

(c)  7f2f*(C)  = 7e(?fL(x(C),u(e),w(€),€)T)  . 

The  problems  in  subsequent  Sections  are  solved  using  the  logarithmic- 

quadratic  loss  penalty  function  W(x,C,r)  defined  as 

P 

W(x,C,r)  t f(x,C)  - r £ In  g,(x,0  + (i/2r)  2 h,(x,C)*,  (1) 

i-1  .5-1  0 

The  following  result  was  obtained  by  Fiacco  (1973)  as  Theorem  3.1. 
TK30RKM  3*  (Approximation  of  first  order  sensitivity  results  and 
determination  of  estimates  from  W(x,C,r)) 

If  assumptions  Al,  A2,  A3  and  A^  hold  for  Problem  P(€),  then  for  (C,r) 
near  (0,0),  there  exists  a locally  unique,  once  continuously  differen- 
tiable  vector  function  y(€,r)  - (x(€,r) ,u(€,r) ,w(€,r))  satisfying 
7xL(x,u,w,£)  - 0, 

ui6i(x,€)  « r,  1*1 

hj(x#0  - w jr, 

with  (x(0,0),u(0,0),w(0,0))  * (x*,u*,w*),  and  such  that  for  any  (C, r) 
near  (0,0)  and  r > 0,  x(f,r)  is  a locally  unique  unconstrained  minimising 
point  of  W(x,£,r),  with  hi(x(€,r),f)  >0,  iHL,,,,,in,  and  Vx2W(x(C ,r) , f ,r) 
is  positive  definite. 

Since  the  system  of  equations  in  Theorem  3 is  identically  equal 
to  zero  at  r * 0,  it  follows  ttet  V^y(f,r)  can  be  calculated  for  (C,r) 
near  (0,0),  The  following  result  was  shown  by  Fiacco  (1973)  following 
his  Theorem  3.1. 

COROLLARY  3.1*  (Convergence  of  estimates  using  W(x,f,r)) 

If  assumptions  At,  A2,  A3  and  A^  hold  for  Problem  1*(C),  then  for  any  < 
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near  0, 
(a) 
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lim  y(€,r)  = y(C,0)  = y(€)f  the  Kuhn-Tucker  triple  characterized 
r-J‘Cf 

in  Theorem  1 ; and 

I'd)  lira  + Vcy(e,r)  = 7fy(€,0)  = V y(C). 

r-K> 

Armacost  and  Fiacco  (1974)  reported  computational  experience 
with  sensitivity  analysis  in  four  sample  nonlinear  programming  problems. 

The  algorithm  uses  the  fact  that  the  Hessian  of  the  penalty  function  is 
positive  definite  for  r small  enough  and  that  the  gradient  of  the  penalty 
function  is  identically  zero  in  a neighborhood  of  the  solution  point.  Thus, 
the  gradient  of  the  solution  point  taken  with  respect  to  the  parameter 
vector  C is  estimated  as 

V x(C,r)  - -V  2W(x,€,r)*17e  2W(x,e,r)|  . (2) 

I x«x(6,r) 

Using  the  fact  that  u^Cpr)  « r/gi(x(€,r) , 6)  and  w^(€,r)  « (l/r)h^(x(Cfr) ,C), 
the  chain  rule  can  be  applied  to  obtain  7^ui(C,r)  and  7^w^.(€fr)  as  shown 
by  Fiacco  ( 1973) • Convergence  was  shown  by  Fiacco  (1973)  following  his 
Corollary  3*1.  The  above  approach  is  equivalent  to  calculating  7^y(€,r) 
directly  from  the  system  of  equations  in  Theorem  3* 

The  logarithmic-quadratic  loss  penalty  function  can  also  be  used 
to  provide  estimates  of  the  first  and  second  order  sensitivity  of  the 
optimal  value  function.  Let  the  optimal  value  penalty  function  be  defined 
as  W*(C,r)  *•  W(x(C,r),C,r),  The  first  order  portion  of  the  sensitivity 
results  developed  by  Armacost  and  Fiacco  (1973)  in  their  Theorem  4 and 
Corollary  4,1  follow, 

THEORFJi  4 * (First  order  sensitivity  of  W*(C#r)  and  estimates  for  f*(0) 

If  assumptions  A1,A2,  A3  and  A4  hold  for  Froblem  P(C)f  then  for  (tpr) 


near  (0,0)  and  r > o, 
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(a)  lim  W*(f,r)  = L*(C)  = f*(C)j 

T-*C 

(b)  V,W*(C,r)  *=  7,L(x,u,w,f)|  (and  (3) 

l(x,utw)=(x(C,r)  ,u(€,r),w(€,r)) 

(c)  lira  +V  W*(f,r)  -V€L(x(C),u(€),w(€),€)  - f»(€). 

r-*0 

Another  estimate  of  the  optimal  value  function  is  obtained  as 

u 

f (C,r)  s f(x(C,r)  ,<:),  Direct  application  of  the  chain  rule  for  differ- 
entiation then  yields  an  estimate  of  the  first  order  sensitivity  of  the 
optimal  value  function  as 

V^e.r)  “ V(x'c)  V(€'r)  + W 

Under  the  given  assumptions,  continuity  assures  that  r(C,r)  -»  f*(C) 

and  V€/(Ctr)  -*-7ef*(€)  as  r 0*.  Thus,  both  V^f^(C#r)  and  V^W*(t',r) 

are  estimates  of  7^f-x(€)  for  r sufficiently  small.  It  is  beyond  the 

scope  of  this  Section  to  explore  the  relationship  between  these  estimates. 

It  is  easily  shown,  however,  that 

u m 

V*(C,r)  = Vef''(t,r)  - £ ^(V^V^C.r)  + 7^) 

P 


+ + ^ho} 


x=x(f,r) 

It  is  easily  shown  that  the  terms  in  the  summations  on  the  right  approach 
zero  as  r approaches  zero,  Armacost  and  Fiacco  (19?4)  used  equations  (2) 
and  (4)  to  examine  the  trajectory  and  convergence  properties  of  the  gradients 
of  the  solution  joint  and  the  optimal  value  function  from  a computational 
point  of  view.  Here,  equation  ('})  is  also  used  to  estimate  the  first 
order  sensitivity  of  tne  optima]  value  function  and  lias  the  advantage 
that  Vrx(f,r)  need  not  l«e  calculated. 
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3,  User  Options  and  Computer  Codes 

The  basic  SUMT-Version  4 computer  program  and  instructions  for 
its  use  are  described  in  My lander,  Holmes  and  McCormick  (1971).  The 
basic  sensitivity  analysis  subroutines,  user  instructions,  and  instructions 
for  integrating  the  sensitivity  package  with  the  SUMT-Version  4 code  are 
described  in  Armacost  and  My lander  (1973).  Briefly,  the  conduct  of  a 
sensitivity  analysis  is  controlled  by  the  variable  NEX0P3  which  is  given 
a value  on  the  "Second  Option  Card"  in  the  SUMT  input  data  deck.  There  are 
four  choicest  no  sensitivity  analysis,  a sensitivity  analysis  at  the  final 
subproblem,  a sensitivity  analysis  at  each  subproblem  along  the  penalty 
function  minimizing  trajectory,  or  a sensitivity  analysis  at  the  final 
subproblem  for  a range  of  differencing  increments.  In  conjunction  with  this 
option,  two  additional  options  are  added  here  and  come  into  play  whenever 
a sensitivity  analysis  is  conducted.  The  first  option  (technically 
Option  4)  is  controlled  by  the  variable  NEX0F4  and  determines  whether  the 
partial  derivatives  of  the  Lagrange  multipliers  will  be  calculated.  When 
the  calculation  is  done,  the  formulas  described  by  Fiacco  (1973)  are  used. 

The  second  option  added  here  (Option  3)  permits  a screening  of  the  parameters 
to  reduce  the  number  of  partial  derivatives  which  are  estimated  by  limiting 
further  analysis  to  those  parameters  which  will  affect  the  optimal  value 
of  the  objective  function  by  an  amount  exceeding  0,1  percent  of  its  current 
value.  This  option  is  controlled  by  the  variable  NEX0P5#  The  estimate  of 
sensitivity  of  the  optimal  value  function  with  respect  to  a particular 
parameter  under  this  option  is  calculated  using  the  Armacost  and  Fiacco 
(1975)  result  which  involves  the  partial  derivative  of  the  Lagranglan  taken 
with  respect  to  the  parameter  under  consideration.  Subroutines  LMULT  and 
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PRESEN  and  related  coding  in  Subroutine  SENS  implement  Option  4 and 
Option  5»  respectively.  Subroutines  SENS,  LMULT  and  PRESEN  are  listed 
in  Appendix  A,  Specific  instructions  for  using  these  two  options  in 
conjunction  with  the  "Second  Option  Card”  are  given  below  in  Table  1, 
This  information  should  be  added  to  Table  5 in  Mylander,  Holmes  and/fVf 
McCormick  (1971). 


TABLE  1 

THE  SECOND  OPTION  CARD 


Option 

Column 

Value 

Meaning 

4 

28 

®0 

Do  not  estimate  the  partial  derivatives 
of  the  estimates  of  the  Lagrange 
multipliers. 

-1 

Estimate  the  partial  derivatives  of  the 
estimates  of  the  Lagrange  multipliers 
whenever  a sensitivity  analysis  of  the 
solution  point  is  conducted. 

5 

35 

-0 

Estimate  the  partial  derivatives  of  the 
optimal  value  function  and  eliminate 
those  parameters  which  do  not  affect  the 
optimal  value  function  from  subsequent 
sensitivity  calculations. 

-1 

Estimate  the  partial  derivatives  of  the 
optimal  value  function  with  respect  to 
all  parameters,  but  continue  all  subse- 
quent sensitivity  calculations  with 
respect  to  all  parameters. 

-2 

Do  not  estimate  the  partial  derivatives 
of  the  optimal  value  function  first. 
Conduct  the  sensitivity  analysis  with 
respect  to  all  parameters. 

A potential  user  of  these  sensitivity  subroutines  should  be  aware 
that  the  penalty  function  coded  in  SUMT-Version  4 does  not  have  the 
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factor  "V"  in  the  quadratic  loss  term  (see  equation  (l)  in  Section  2), 
Therefore,  the  expressions  for  the  several  gradients  have  an  additional 
factor  of  "2’'  in  the  comjjuter  program  which  does  not  appear  in  the 
supporting  theory  of  Section  2. 


4,  New  Computational  Experience 

In  this  Section,  the  four  sample  problems  of  Armacost  and  Fiacco 
(1974)  are  examined.  Specifically,  the  convergence  of  the  partial  deriva- 
tives of  the  Lagrange  multipliers  is  examined  and  the  estimates  of  the 
gradient  of  the  optimal  value  function  obtained  by  the  chain  rule  (equation 
(4))  and  by  the  gradient  of  the  Lagrangian  (equation  (3))  are  compared. 

The  problems  are  designated  by  the  same  letters  as  in  the  original  paper. 
Consider  first  a simple  convex  program.  The  problem  is 

minimize  f(x,C)  = x^  + 

2 2 2 ® 
subject  to  g^(x,C)  « - x^  >0, 

for  > 0.  The  analytical  solution  point  and  its  gradient  are  given  in 

Armacost  and  Fiacco  (1974)  as 


x(€) 


f*(0 


v(<0 


i + e. 
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(1  + V) 


(1  ♦ e2  ) 


V€f*(€)  » (af»(6)/861(  8f*(€)/d€  ) 


“ V 1 + €2  , -€j€2/  VT  + €22^  , 


The  Lagrange  multiplier  and  Its  gradient  are  analytically  determined 
to  be  

. 4 / 4 X /-  2 / n/. 


u*(€)  - yi^2  /2€t  , 


7€u*(€)  - (-V1  + €22  /2Z*9  6a/(2€tV  1 * ^2)  ). 

The  numerical  example  had  - 2 and  €2  - 1 yielding  the  following 
numerical  results t 


f*  - -2V2  , 

V€f*  - 

(-V2.  -VI ) , 

-V t 

-V2/2  V2/2] 

**  “ , 
1-V2J 

V^X*  - 

.-V2/2  -Vl/d’ 

u*  - V2A 

V€u*  - 

(-V2/8.  \^/8) 

« 0.353  , 

(-.177,  .177)  . 

The  numerical  results  obtained  by  the  computer  program  are  In* 
eluded  In  Table  2 for  the  optimal  value  function  and  Lagrange  multiplier 
sensitivity.  The  values  of  the  first  order  optimal  value  function  sensi- 
tivity computed  both  by  the  chain  rule  (equation  (4))  and  by  taking  partial 
derivatives  of  the  Lagranglan  with  respect  to  the  parameters  (equation  (3)) 
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TABLE  2 

TRAJECTORY  RESULTS  FOR  PROBLEM  B 


Lagrangian  Chain  rule 

Subproblem  f df/de.,  of/o^ 

u 

bu/d€1 

4u/4€2 

1 -1.9999  -1.9999  - .9999  -1.3333  -1.3333 

.4999 

-.3333 

.1666 

2 -2.5393  -1.5^0  -1,2947  -1.4087  -1.4088 

.3860 

-.2100 

.1761 

3 -2.7765  -1.4439  -1.3833  -1.4139  -1.4139 

.3609 

-.1844 

.1767 

4 -2.8128  -1.4243  -1.4064  -1.4142  -1.4142 

.3560 

-.1790 

.1768 

5 -2.8245  -1.4127  -1.4123  -1.4142  -1,4142 

.3531 

-.1768 

.1768 

6 -2.8274  -1.4137  -1.4137  -1.4142  -1.4142 

.3532 

-.1767 

.1768 

7 -2.8282  -1.3899  -1.4141  -1.4142  -1.4142 

.3475 

-.1737 

.1768 

Analytical  -2.8282  -1.4142  -1.4142  -1.4142  -1.4142 

.3537 

-.1769 

.1768 

are  presented  In  parallel.  The  results  are  also  plotted  in  Figure  i and 
portray  the  type  of  convergence  experienced.  While  the  previous  results 
toy  Armacost  and  Fiacco  (1974)  clearly  indicated  a stability  of  the  solution 
point  and  optimal  value  function  and  their  gradients  taken  with  respect  to 
the  parameters,  Table  2 indicates  that  with  the  Lagrange  multipliers,  the 
same  sort  of  stability  is  not  found.  It  is  well  known  that  with  barrier 
functions,  the  estimates  of  the  Lagrange  multipliers  decrease  in  accuracy 
as  the  boundary  is  approached.  The  change  in  the  value  of  u between  sub- 
problems  6 and  7 is  an  Indication  of  this.  It  is  no  surprise,  therefore, 
that  the  estimates  of  the  gradient  of  the  Lagrange  multiplier  behave  in 
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a similar  way.  In  addition,  since  the  estimate  of  the  gradient  of  the 
optimal  value  function  obtained  by  evaluating  the  partial  derivatives  of 
the  Lagranglan  taken  with  respect  to  the  parameters  Includes  the  estimates 
of  the  Lagrange  multipliers,  it  too  will  not  be  as  accurate  an  estimate 
as  the  boundary  is  approached. 


The  next  problem  considered 

is 

a nonconvex  program  with 

constraint.  The  problem  is  to 

minimize  f(x,€) 

- 

X1  * x2  * 

In  Xj  - x^ 

subject  to  g^(x,€) 

m 

V * *2 

> 0, 

«2(x.€) 

- 

X1 

> of 

«3(m) 

m 

x3  ~€1 

> 0, 

ht(x,e) 

- 

2 2 
x3  * x4 

- S*  - «• 

where  ^ > 0 and  €2  > 0, 

The  analytical  solution  isi 

f*(€)  - In  -W 

- c 

2 

1 

U*(€)  - ^7 

x1*(f)  ; x,«(C)  - 0, 

Y < 

3 

(0  ii  * 

Uj*(0  - u2*(0  * 1, 

u3*(€)  -l/€t  ♦ €t/  V €22  - t*  , 

and  «1*(€)-i/(2V  7^  - 7^  . 

The  numerical  example  used  has  - 1 and  ^ • 2 . The  numerical  solution 
derived  analytically  for  the  solution  point,  Lagrange  multipliers  and  their 
gradients  area 

V*  * (1.578,  -1.157)  , 


f*  * -1.732 


x*  a 


' 0 " 

— • 
0 

0 

0 

1 

*5 

a 

0 0 
1 0 

.1.732, 

r.577  1.154, 
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* « 
1 

0 

0 

u*  w 

1 

.1.732 

. V€u*  * 

0 

-.231 

0 

-.38 5. 

v*  « .289  , V€w*  a (.0962,  -.1925)  . 

The  computed  results  for  a trajectory  sensitivity  analysis  are  shown  in 

•a 

Tfcble  3. 

The  partial  derivatives  of  the  Lagrange  Multipliers  shown  in  Table  3 
are  the  only  ones  which  are  non-zero.  Note  that  at  the  fifth  and  sixth 
subproblems,  the  estimates  of  the  gradients  of  u^  and  w^  are  reasonably 
close  to  the  true  values  determined  analytically.  It  is  at  that  point  also 
that  the  estimates  of  u^  and  w^  are  the  closest  to  their  true  values  at 
the  solution  point.  Notice  also  that  the  estimate  of  df*/^  obtained 
from  the  partial  derivative  of  the  Lagranglan  with  respect  to  is  reason- 
ably close  to  the  true  value.  It  is  entirely  dependent  on  w^  and  as  the 
estimate  for  w^  becomes  less  accurate , the  error  will  be  reflected  in 
the  first  order  sensitivity  estimate  of  the  optimal  value  function  taken 
with  respect  to  In  the  following  example,  the  need  for  careful  attention 
to  the  differencing  increment  is  illustrated  when  the  parameters  are  the 
right-hand  sides  of  the  constraints.  For  Problem  C,  a sensitivity  analysis 
was  conducted  at  the  final  subproblem  for  a range  of  differencing  increments. 

The  results  were  that  the  sensitivity  estimates  remained  fairly  constant 

-?  -12 

over  the  range  of  differencing  increments  from  10  to  10  , Thus,  the 

source  of  error  in  this  case  is  solely  the  lack  of  accuracy  of  the  estimates 
of  the  Lagrange  multipliers  for  the  binding  constraints. 

Two  related  problems  called  the  Shell  Primal  and  the  Shell  Rial 
were  presented  by  Anaeoet  and  Fiaeoo  (19?**).  However,  computational 
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results  were  presented  only  for  the  Shell  Dual,  (The  Shell  Dual  was 
developed  as  a test  problem  by  the  Shell  Development  Company  and  used  by 
Colville  (1968)  in  his  comparative  analysis  of  nonlinear  programming 
codes.)  Computational  results  are  presented  below  for  both  the  Shell 
Primal  and  the  Shell  Dual.  The  first  problem  considered  is  the  Shell 
Primal. 

n n n n ~ 

minimize  f(x,€)  - £ e,x.  ♦ Z Ex,c.,x.  t Ed.x/* 

j-1  3 3 i-1  j-1  1 13  3 j-1  3 3 


subject 


H 

to  ^(x,*)  - Z ai jXj  - >0,  i-1,..., 

J— 1 


with  Xj  > 0,  j-l,.,.fn.  The  dual  problem  is  much  more  difficult  to  solve 
and  is  the  one  most  often  used  In  computational  comparisons.  The  Shell 
foal  is  to 


n n 


maximize  f(x,€)  • Z €jr . - Z Z x.c  .x.  - 2 Z d.x  - 

J-1  x 3 i-1  j-1  1 1 3 3 i-1  1 1 


subject  to  g1(x,€)  - ♦ 2 Z c^x^  ♦ 3*^^ 

m 

* Z a , ,y . > 0,  1— l,».4»ny 
j-1  31  3 


with  x1  > 0,  1-1,. ..,n,  and  y^  > 0,  j-1, In  the  numerical  example 
used  here,  n - 5 and  m - 10.  The  problem  data  is  given  in  Table  3 of 
Armacoet  and  Flacco  (197*0  and  in  Appendix  C.  The  parameters  of  the 
sensitivity  analysis  are  the  variables  1-1,,,,, 10,  the  components  of 
the  right-harvi  sides  of  the  primal  constraints. 

As  a urlef  aside,  the  computational  solution  of  the  Shell 
.'Mtl  provided  the  motivation  for  some  ^f  the  recent  work  in  parametric 
sensitivity  analysis  by  Armacost  nd  Fiacco.  Specifically,  Armacost  and 
Klaccn  (19'A)  noted  that  in  solvin*,  vhe  dual  problem,  the  partial 
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derivatives  of  the  dual  variables  with  respect  to  the  right-hand  sides 
of  the  primal  constraints  were  obtained.  With  the  correspondence  between 
the  dual  variables  and  Lagrange  multipliers  and  their  interpretation  as 
the  partial  derivatives  of  the  optimal  value  function  with  respect  to 
the  right-hand  sides  of  the  primal  constraints,  it  apptared  that  the 
second  order  partial  derivatives  of  the  optimal  value  function  had  been 
obtained.  The  calculations  supported  this  conjecture  since  the  matrix  of 
partial  derivatives  of  the  dual  variables  with  respect  to  the  parameters 
was  symmetric.  Arina  cost  and  Fiacco  (1976)  have  shown  that  when  the 
parameters  are  the  components  of  the  right-hand  side  only,  then  the 
gradient  of  the  Lagrange  multiplier  vector  taken  with  respect  to  the 
parameters  is  the  Hessian  of  the  optimal  value  function.  This  matrix  will 
be  computed  using  the  Shell  Primal  and  then  compared  with  the  Hessian 
obtained  by  solving  the  Shell  Dual, 

In  solving  all  of  the  sample  problems,  the  option  to  screen  the 
sensitivity  estimates  was  used  resulting  in  the  partial  derivatives  being 
computed  only  for  those  parameters  which  affected  the  optimal  value 
function  by  more  than  O.i  percent  of  its  current  value.  Annotated  com- 
puter output  for  the  final  subproblem  with  sensitivity  analysis  data  for 
the  Shell  Primal  is  shown  in  Figure  2,  (The  annotation  applies  to  the  com- 
puter output  in  Figures  3#  5 and  ' as  well,)  The  parameters  are  represented 

by  the  letter  *’AM  vice  In  the  computer  output.  Similar  output  for  the 
Shell  Dual  ir.  siiown  in  Figure  -,  ( Compare  the  sensitivity  analysis  portion 
with  -Lure  **  of  Armaco^t  and  Fiacco  (i >?**).)  The  sensitivity  estimates  for 
both  problems  were  obtained  by  conducting  a trajectory  sensitivity  analysis, 
i.e.,  a sensitivity  analysis  performed  at  each  subproblem  along  the 
minimising  trajectory,  Cince  the  Ci *11  :\isl  is  a maximisation  uroblem 
and  51‘r'T  is  coded  to  solve  a minimi  sat  ion  problem,  Problem  K is  solved 
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-Shell  Prlcal  subproblen  and  sensitivity  output 
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\ 

ij 

, 


Annotation 


Identifier 

1 

2 

3 


4 


5 


6 


7 


8 


Ci 

y' 


10 

11 

1? 


13 


Meaning 

F = the  value  of  the  objective  function 
at  the  current  solution  point, 

The  value  of  the  components  of  the  solution 
point  of  the  current  subproblem,  here, 

x^ , , , * , x^. 

The  value  of  the  constraints  evaluated  at  the 
current  solution  point,  i.e,, 

The  data  corresponding  to  1 - 3 above  when 
the  solution  point  is  a second  order  estimate 
based  on  the  values  in  2, 

The  same  as  4,  but  the  solution  point  is  a 
first  order  estimate  based  on  the  values  in 
2.  These  values  and  those  in  4 are  extrapo- 
lations from  the  current  solution  estimate. 

The  estimates  of  the  Lagrange  multipliers  based 
on  the  current  value  of  r (RHO)  and  the  current 
estimates  of  the  solution  point,  in  this  case 
as  ui  « r/gi(x,€),  i^l,.,.,lO. 

The  value  of  the  parameters,  C^,,  i*=l,,.,,10. 

The  value  of  the  differencing  increment  used 
in  the  central  differencing  formula  for  each 
of  the  parameters. 

The  estimates  of  the  gradient  of  the  optimal 
value  function  calculated  by  equation  (3). 

The  parameters  whose  associated  partial  derivatives 
affect  the  optimal  value  function  by  more  than 
0,001  of  its  current  value. 

The  first  order  sensitivity  of  the  solution 
point  calculated  by  equation  (2), 

The  first  order  sensitivity  of  the  Lagrange 
multipliers  calculated  using  the  method 
described  following  equation  (2), 

The  partial  derivative  of  the  optimal  value 
function  taken  with  respect  to  the  indicated 
parameter  and  calculated  ly  equation  (4), 


Fir.,  2. — Continued, 
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Pig#  3# — Shell  Dual  subproblem  and  sensitivity  output. 
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by  minimizing  the  negative  of  the  objective  function  and  therefore,  the 
results  in  Fig*  3 for  the  objective  function  value  and  the  values  of 
its  partial  derivatives  must  be  multiplied  by  -1  to  obtain  the  correct 
results.  The  variables  X(l)  - X(10)  in  the  Shell  Dual  correspond  to  the 
Lagrange  multipliers  in  the  Shell  Primal  and  the  variables  X(ll)  - X(15) 
in  the  Dual  correspond  to  X(l)  - X(5)  in  the  Primal,  The  components  of 
the  Hessian  of  the  optimal  value  function  then  are  the  partial  derivatives 
of  X(l)  - X(10)  in  the  Shell  Dual  and  the  u-derivatives  in  the  Shell  Primal, 
The  dual  variables,  Lagrange  multipliers,  the  partial  derivatives  of  the 
optimal  value  function  obtained  by  means  of  the  Lagrangian,  and  those 
obtained  using  the  chain  rule  are  compared  in  Table  4, 

TABLE  4 

FIRST  ORDER  SENSITIVITY  COMPARISON 
Shell  Primal  Shell  Dual 

i af/ae1(Lag.)  af/ae^cX)  af/ae^iag.)  af/ae^c.R.) 


1 

,l68xl0-5 

.I68xl0~5 

- 

,168xl0"5 

1 

0.0 

- 

2 

.312x10 

.312x10 

- 

.313x10“* 

0.0 

- 

3 

5.1649 

5.1649 

5.1740 

5.1742 

5.1741 

5.1741 

4 

.437X10"4 

.437X10"4 

- 

.438X10*4 

0.0 

- 

5 

3.0555 

3.0554 

3.0610 

3.0610 

3.0610 

3.0612 

6 

11.8191 

11.8190 

11.8395 

11,8405 

11.8406 

11.8397 

7 

,159xl0“5 

.159x10'“’ 

- 

,159xl0'5 

0.0 

- 

8 

.107xl0*5 

,107xl0'5 

- 

.10?xl0'5 

0.0 

- 

9 

.1046 

.1046 

.1038 

.1043 

.1044 

.1040 

10 

.889x10"^ 

.809X10-4 

- 

,891xl0J* 

.177xl0'3 

Since  the  preliminary  screening  option  was  used,  detailed  sensitivity 
estimates  for  the  parameters  which  correspond  to  the  non-binding  Inequality 
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constraints  were  not  computed  in  Figures  2 and  3*  In  the  Shell  Primal, 
the  Lagrange  multipliers  and  the  partial  derivatives  obtained  by  means  of 
the  Lagrangian  correspond  exactly.  As  noted  in  the  discussions  following 
Problems  B and  C,  the  estimates  of  the  Lagrange  multipliers  for  the 
binding  constraints  are  very  sensitive  near  the  boundary  of  the  constraint 
set.  This  explains  the  slight  variation  between  these  estimates  and  the 
other  sensitivity  estimates  shown  in  Table  4 for  the  binding  Primal 
constraints. 

Now  compare  the  second  order  sensitivity  estimates  with  respect 
to  the  parameters  which  are  the  right-hand  sides  of  the  binding  Primal 
constraints.  Let  and  Hp  denote  the  submatrices  of  the  Hessian  of  the 
optimal  value  function  for  the  Dual  and  Primal  respectively,  obtained  by 
deleting  the  rows  and  columns  corresponding  to  the  non-binding  Primal 
Inequality  constraints.  Thus,  the  components  of  Hp  are  dx^dCj,  i*  j“3»5*6»9, 
and  the  components  of  Hp  are  du i» j“3*5f6,9.  Fro®  the  computer 
output,  these  matrices  arei 


4,0642 

-.5653 

3.4650 

.1894 

-.5653 

.5043 

-.6151 

.0237 

3.4650 

-.6151 

7.1850 

1.5885 

. .1894 

.0237 

1.5885 

.8343. 

'4.0710 

-.5663 

3.4715 

.1899 

-.5662 

.5051 

-.6158 

.0239 

3.4715 

-.6158 

7.1929 

1.5890 

. .1899 

.0239 

1.5890 

.8347. 

Both  Hp  and  Hp  are  symmetric,  and  while  the  agreement  is  not  exact,  it  is 
very  close  as  anticipated, 

Arm&cost  and  Fiacco  (1974)  indicate  that  the  differencing  interval 
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used  in  the  sensitivity  analysis  can  affect  the  accuracy  of  the  results. 

For  this  reason,  an  option  was  provided  by  Armacost  and  Mylander  (1973) 
to  conduct  the  sensitivity  analysis  at  the  final  subproblem  for  a range 
of  differencing  intervals.  It  is  even  more  important  to  be  cautious 
when  dealing  with  right-hand  side  perturbations  as  the  following  discussion 
indicates.  The  results  shown  above  for  the  Shell  Primal  and  the  Shell 
Dual  were  obtained  with  a trajectory  sensitivity  analysis  and  at  the 
final  subproblem,  the  differencing  increment  used  in  the  central  differ- 
encing formulas  was  lo”1*.  The  Hessian  subna trices  Hp  and  Hp  were  found 
to  be  very  close.  When  the  problems  were  solved  using  a sensitivity 

analysis  at  the  final  subproblem  only  with  a differencing  increment  of 
-9 

10  , the  diagonal  elements  of  Hp  were  considerably  different  from 

those  of  Hp,  The  problems  were  solved  again  with  a sensitivity  analysis 
performed  at  the  final  subproblem  for  a range  of  values  of  the  differencing 
increment,  ranging  from  10“^  to  10"*1.  The  components  of  Hp  remained 
constant  as  did  the  non-diagonal  elements  of  Hp  which  were  equal  to  the 
non-diagonal  elements  of  Hp.  The  diagonal  elements  of  Hp  did  not  remain 
constant  and  their  variation  is  depicted  in  Table  5* 

The  final  example  considered  in  this  section  is  called  the  cattle 
feed  problem.  It  was  formulated  and  originally  presented  by  van  de  Panne 
and  Popp  (1963)*  Armacost  and  Fiacco  (197*0  presented  the  cattle  feed 
problem  to  illustrate  an  application  of  the  sensitivity  analysis.  The 
additional  sensitivity  results  are  presented  here  for  completeness.  The 
problem  is  a chance-constrained  program  to  determine  the  mix  of  inputs  to 
cattle  feed  that  will  satisfy  nutritive  constraints  and  minimise  the  cost 
of  the  cattle  feed.  The  protein  content  of  the  components  is  a random 
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TABLE  5 


VARIATION  IN  THE  COMPONENTS  OF  Hp 


A 

&2f/df 

2 

1 

i-3 

i*5 

i“6 

1-9 

io4 

-3148.3 

-383.81 

-89157.9 

.8342 

1(T7 

-27.228 

-3.328 

-851.345 

.8347 

to'8 

3.758 

.466 

-1.389 

.8347 

1<T9 

4.067 

.504 

7.107 

.8347 

o 

o 

4.071 

.505 

7.191 

.8347 

ID"11 

4.071 

.5 05 

7.192 

.8347 

variable,  normally  distributed  with  a mean  and  variance  determined  experi- 
mentally. The  application  of  the  sensitivity  analysis  included  the 
standard  deviations  as  parameters  with  the  interpretation  that  if  the 
solution  were  sensitive  to  a standard  deviation,  more  sampling  would  be 
indicated  in  order  to  obtain  a sharper  estimate.  The  statement  of  the 
problem  is 


minimize 
subject  to 


f(x,€)  - CjXj  + CjXg  + CyXj  + ckxk 

gjtx.e)  - a1*1  + EjjXj  + aycj  + a4xu  - ^ > 0, 
k2(x.€)  - ^ ♦ MjXj  + Mux4  F 


hjtx.O  - \ + Xj  x3  + x^  - 1 - 0, 

with  x^  > 0,  i-1,2,3,**.  The  notation  is  slightly  different  from  Armacoet 
and  Fiacco  (197**).  Here,  the  parameters  in  the  sensitivity  analysis  are 
denoted  The  correspondence  with  the  previous  work  is  6^  • O^,  i-1,2,3,1*, 
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“ v and  Cr,  - o , The  problem  data  from  Table  6 of  Armacost  and 

5 rt  ' 1 m / m r 

Fiacco  (ly?4)  are  included  in  Appendix  C,  The  sensitivity  results  for 
Problem  F are  shown  in  Figure  Again,  the  preliminary  screening  option  wa 
used  to  avoid  calculating  sensitivity  estimates  which  did  not  affect  the 
optimal  value  function  significantly.  Here,  the  components  of  the  gradient 
of  the  optimal  value  function  obtained  directly  from  the  gradient  of  the 
Lagrangian  and  the  estimate  of  the  gradient  obtained  by  application  of  the 
chain  rule  are  very  close.  The  values  of  the  Lagrange  multipliers  estim- 
ated in  the  SUHT  program  are  u^  «=  0.58037,  u^  **  -.41005  and  w^  = -18,3738, 
The  multipliers  and  u 2 correspond  to  and  df*/^  respectively. 

Also  note  that  second  order  sensitivity  information  for  and  is 
available  since  they  are  the  right-hand  sides  of  the  two  inequality  con- 
straints, Namely, 

a2f*/d€?2  - ou1/i€?  - .0112929, 

32f*M6a€7  " Sui/ae6  “ •45?06x10-6  , 

32f*/3€?3€6  - 3u2/3€?  - .45671xl0’6  , 
and  32f*/3€62  = '^'ai6  “ • 4621 2x1 °'3  . 

The  results  shovm  in  Figure  4 were  obtained  using  a sensitivity  analysis 
at  the  final  subproblem  for  a range  of  differencing  increments.  For 
other  values  of  the  differencing  increments  not  shown  here,  the  optimal 
value  function  sensitivity  estimates  with  respect  to  6^  and  obtained 
directly  from  the  gradient  of  the  Lagrangian  and  the  cross-port ials  of 
u ^ and  with  respect  to  and  do  ,fary  somewhat  while  the  other 
sensitivity  estimates  remain  relatively  constant.  In  a practical  sense, 
however,  the  cross-part  ials  of  u ^ and  u2  with  respect  to  and  are 
constant  since  they  are  of  the  order  of  10*^. 
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5«  Large-scale.  Multi- Item  Inventory  Model 

Traditionally,  Inventory  models  have  been  formulated  to  minimize 
some  function  of  the  ordering,  holding  and  shortage  (or  backorder)  costs 
subject  to  various  constraints,  Schrady  and  Choe  (1971)  have  formulated  an 
inventory  model  which  appears  to  have  much  greater  relevance  for  an  inventory 
system  such  as  the  U,  S.  Naval  supply  system.  For  the  Navy,  the  costs  used 
In  the  traditional  models  may  be  quite  artificial  while  the  real  objective  of 
the  system  is  to  maximize  the  service  to  the  Fleet,  an  objective  equivalent 
to  minimization  of  stockouts.  In  addition,  the  stock  points  of  the  Naval 
supply  system  have  investment  and  reorder  workload  constraints  that  are  real 
and  binding,  Schrady  and  Choe  consider  a multi-item  inventory  system  with  the 
specific  objective  function  of  minimizing,  the  total  time-weighted  shortages. 
The  decision  variables  are  the  "reorder  quantities"  and  the  "reorder  points," 
the  decisions  of  how  much  to  order  and  when  to  order,  for  each  item  in  the 

inventory.  Clearly,  it  is  of  interest  to  know  "if"  and  ":y  how  much"  these 
variables  and  the  value  of  the  objective  function  will  change  if. certain 
parameters  change,  Schrady  and  Choe  solved  a small  example  problem  using 
the  SUMT  computer  code,  KcCormick  (1972)  has  shown  how  the  special  structure 
of  this  inventory  model  can  be  used  to  facilitate  the  use  of  the  SUHT  code  to 
solve  very  large  inventory  problems.  Here,  the  sensitivity  analysis  is 
applied  to  this  model,  a large-scale  inventory  system,  and  it  is  illustrated 
by  means  of  the  example  used  by  Schrady  and  Choe, 

The  model  presented  here  is  due  to  Schrady  and  Choe,  An  extension  of 
the  model  by  KcCormick  (1972)  includes  constraints  on  storage  volume  and 
probability  of  depletion.  The  sensitivity  results  are  easily  applied  to  the 
extended  model.  Several  aaeumptlons  specify  the  nature  of  the  model.  The 
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first  is  that  all  demand  which  occurs  when  the  on-hand  stock  is  ssrof  is  back* 
ordered*  The  model  is  probabilistic  in  that  the  lead  tine  demand  is  a random 
variable.  Specifically,  for  the  1th  variable,  it  is  assumed  that  the  demand 
which  occurs  during  the  time  between  the  placement  of  an  order  and  its  receipt 
ty  the  stock  point  is  normally  distributed  with  mean  and  variance  , 

For  th,  1th  Item,  let 

c^  - item  unit  cost  (in  dollars), 

- mean  demand  per  unit  time  (in  units), 
r^  - reorder  point, 

4^  - reorder  quantity, 

<|>(x)  - the  Konaal(0,l)  density  function,  and 

Q(s)  * the  Normal ( 0,1)  complementary  cumulative 
distribution  function  - £V*>  a* . 

In  addition,  let  be  the  investment  limit  in  dollars,  let  1^  be  the 
reorder  workload  limit  and  let  N be  the  total  number  of  items  in  the  inventory. 
The  detailed  development  of  the  model  is  omitted  here  and  the  reader  is  referred 
to  Schrady  and  Choe  (1971)  or  NcCoxmlck  (1972),  A function  needed  from  that 
development  for  the  final  model  is 

*!<*!>  - Kot2  ♦ (rt  - M1)2)<*K(r1  - 

Then  the  general  multi-item  model  of  Schrady  and  Choe  is 

If 

minimi**  Z(Q,r)  - 
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N 

subject  to  gt(Q,r)  - Kj  - £ + ^/2  " Mj)  > °. 

N 

^(Q.r)  - K2  - I Jj/Jj  > 0 , 

with  unrestricted,  Qj  > 0,  1-1,... ,N,  Q - (Qj,...,Q(,)T  and  r - (r^...,^)1. 

To  jut  the  problem  In  (x,<)  notation,  sake  the  following  Identifications 
for  1-1,...,  Ni 

*21-1  " Ql' 

*21  * rl> 

€m-i  ■ \ 

*41  ‘ °1» 

*41+1  “ cl» 

*41+2  Al* 

and  - K^f  - K2.  Rewrite  Problem  S-C  u 

N 

minimise  f(x,€)  - 2 &i(*2i*€)/*2i-i 

1-1  S-C(€) 

N 

subject  to  SjCx,*)  “ €j  - * *21-V2  ’ *4l-l> 


K2(*.€)  *€2  “ Z* 41+2^*21  *1  - °’ 

>0,  Xj^  unrestricted,  1-1,. ..,»,  and 

^i^i**^  “ K(eui2  * (*21  ’ W^^l  ’ €4l-i^4i) 


“ €Ul/x2i  "*4l-l^  ^^*21  ” *4i«l^*4l^* 

Schrady  and  Choe  consider  a three  Item  sample.  The  problem  data 

mad  the  Initial  starting  point  for  the  SUHT  program  are  shorn  In  Table  6 , 
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TABLE  6 

MULTI-ITEM  INVENTORY  PROBLEM  DATA 


Item 

i«l 

i«2 

1-3 

Data* 

100 

200 

300 

V^I 

100 

100 

200 

ci‘W 

1 

10 

20 

vw 

1,000 

1,500 

2,000 

Starting  point* 

V*21-l 

600 

270 

300 

V’zi 

200 

260 

400 

In  addition,  Kx  » ^ - $8,000  and  Kg  * e,  - 15. 

Figure  5 contains  the  computer  output  for  the  final  subproblem  of  a 
trajectory  sensitivity  analysis  for  Problem  S-C(€)  using  the  data  of  Table  6. 
The  results  Indicate  that  the  optimal  value  function  is  sensitive  to  parameters 
2,  5,  8,  9,  12,  and  13,  i.e.,  Kg,  c1#  Og,  Cg,  Oj  and  respectively.  The 
fact  that  the  solution  is  sensitive  to  the  values  of  the  standard  deviations 
of  the  Dead  time  demand  of  two  items  lets  the  decision  maker  know,  that  since, 
these  parameters  were  obtained  by  sampling,  a possible  action  may  be  to  con- 
duct additional  sampling  in  order  to  sharpen  the  estimate  of  the  standard 
deviation. 

The  solution  value  is  also  very  sensitive  to  all  of  the  item  costs. 

If  the  structure  of  Problem  S-C  is  examined,  this  result  is  most  surprising 
since  the  c^  appear  only  in  the  investment  constraint  and  the  optimal  value 
function  is  not  very  sensitive  to  the  investment  limit  K^,  Recall,  however, 
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fig.  5.  —Computer  output  for  Schrady-Choe  inventory  problem 
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that  when  the  partial  derivatives  are  used  as  sensitivity  estimates,  one  is 
effectively  saying  "how  much  will  the  objective  function  (or  solution  point) 
change  if  the  parameter  is  increased  by  one  unit?"  Suppose  c^  (parameter  5) 
is  increased  by  one  unit  from  one  to  two.  Using  a linear  estimate,  the 
objective  function  is  expected  to  increase  by  2,16  to  15*15  and  (variable 
x^)  is  expected  to  decrease  by  208  units  to  325.  Figure  6 contains  the 
computer  output  for  the  perturbed  problem  with  c^  - 2,  Note  that  the  value 
of  the  optimal  value  function  increased  to  14,99  and  decreased  to  411,  The 
changes  were  quite  large  and  in  the  directions  expected  but  of  course  not 
as  large  as  the  linear  estimate.  Notice,  however,  that  this  change  in  c^ 
represented  a 100  % increase  in  the  value  of  the  parameter.  For  purposes 
of  comparison,  consider  a similar  change  in  the  investment  constraint 
(parameter  1)  to  which  the  solution  is  apparently  insensitive.  Using  the 
optimal  value  function  sensitivity  estimate  (-.00517),  a 1C0£  increase  in 
the  value  of  the  parameter  is  8,000  and  consequently,  a linear  estimate  of 
the  expected  change  in  the  optimal  value  function  would  be  (8,000  x -,00517  “) 
-41,36,  The  above  example  does  not  imply  ttiat  equivalent  percentage  changes 
in  parameter  values  is  relevant  'out  rather  is  meant  to  emphasize  that  the 
sensitivity  estimates  are  valid  for  a small  neighborhood  of  the  given 
parameter  values  and  as  such  represent  instantaneous  changes. 

The  solution  values  in  Figure  5 are  slightly  different  from  those 
presented  in  Schrady  and  Choe  (1971).  This  difference  is  due  to  the  use  of 
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two  different  approximations  to  the  complemehtary  cumulative  normal  distribu- 
tion. While  the  solution  values  vary  only  slightly  with  these  different 
estimates  of  the  normal  distribution,  the  variations  in  the  Lagrange  multipliers 
and  some  sensitivity  estimates  are  greater.  For  the  system  library  subroutine 
using  the  normal  distribution  error  function  (erf)  to  estimate  the  complementary 
cumulative  distribution,  the  minimizing  trajectory  of  subproblems  showed  that 
the  estimate  of  the  Lagrange  multipliers  deviated  from  what  appeared  to  be  a 
relatively  constant  value  as  the  final  subproblmm  was  approached.  As  expected, 
the  optimal  value  function  sensitivity  estimates  obtained  directly  from  the 
Lagrangian  varied  in  almost  direct  proportion  with  the  estimates  of 
the  Lagrange  multipliers.  In  Figures  5 and  6,  the  estimates  of  the  optimal 
value  function  sensitivity  obtained  by  the  two  different  methods  (chain  rule 
and  Lagrangian)  are  in  close  agreement.  This  was  not  the  case  using  the  erf- 
related  subroutine  used  by  Schrady  and  Choe.  This  is  the  same  effect  experienced 
in  other  problems  as  discussed  in  Section  4.  In  this  case,  however,  the 
x-derivatives  are  affected  to  a slight  degree.  The  major  source  of  this 
error  appears  to  be  the  lack  of  necessary  precision  associated  with  the  erf- 
related  subroutine  for  the  complementary  cumulative  normal  distribution.  The 
conclusion  can  only  be  that  one  must  proceed  with  caution.  From  this  and  other 
examples,  however,  it  appears  that  convergence  of  the  Lagrange  multiplier 
estimates  along  the  minimizing  trajectory  is  a good  indication  that  the  sensi- 
tivity estimates  will  be  accurate. 

This  large-scale  inventory  example  illustrates  a potential  Hreal  world" 
application  of  sensitivity  analysis.  It  also  highlights  the  need  for  a 
careful  interpretation  of  the  sensitivity  information  as  well  as  the  recurrent 
call  for  caution  in  the  use  of  a numerical  algorithm. 
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6,  Conclusions 

The  purpose  of  this  paper  was  to  present  examples  of  computational 
implementation  of  sensitivity  analysis  with  respect  to  the  optimal  value 
function  and  the  Lagrange  multipliers  using  the  theoretical  results  of 
Fiacco  (1973)  and  Armacost  and  Fiacco  ( 1975*  1976).  All  of  the  problems 
presented  led  to  new  insights  into  the  computational  aspects  of  this  type 
of  sensitivity  analysis. 

The  major  conclusion  is  that  the  sensitivity  estimates  of  the 
optimal  value  function  obtained  directly  from  the  gradient  of  the  Lagrangian, 
the  partial  derivatives  of  the  Lagrange  multipliers  associated  with  the 
binding  constraints,  and  to  a lesser  extent,  the  partial  derivatives  of 
the  solution  point,  are  dependent  on  the  accuracy  of  the  estimate  of  the 
Lagrange  multipliers  calculated  in  the  penalty  function  algorithm  and 
subsequently  used  in  the  gradient  of  the  Lagrangian.  If  the  estimates  of 
the  Lagrange  multipliers  along  the  minimizing  trajectory  converge  to  a 
common  value,  then  it  appears  that  the  sensitivity  estimates  will  converge 
to  their  true  values  provided  the  differencing  increment  is  satisfactory. 

If,  however,  the  Lagrange  multipliers  converge  and  then  vary,  the  sensitivity 
estimates  should  be  viewed  with  caution.  (Note  that  the  Lagrange  multiplier 
estimates  are  available  in  the  standard  SUMT  output  and  a trajectory  analysis 
is  not  required.) 

Armacost  and  Fiacco  (197*0  noted  that  the  differencing  increment 
used  in  the  central  differencing  formulas  was  a potential  source  of  error 
depending  on  the  scaling  of  the  problem.  The  same  caution  applies  when 
dealing  with  Lagrange  multiplier  and  optimal  value  function  sensitivity. 
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The  example  problems  which  included  right-hand  side  parameters 
indicated  that  the  optimal  value  function  second  order  sensitivity  estimates 
are  themselves  very  sensitive  to  the  Lagrange  multiplier  estimates,  the 
differencing  increment  and  the  scaling  of  the  problem.  It  appears  that 
further  analysis  is  needed  before  this  particular  program  can  be  used  to 
take  advantage  of  second  order  sensitivity  estimates  to  improve  algorithm 
performance.  This  is  particularly  true  for  the  second  partial  derivatives 
of  the  optimal  value  function  taken  with  respect  to  the  right-hand  side 
of  a binding  constraint. 

Computer  time  was  provided  by  The  George  Washington  University 
Computer  Center. 


- M - 


REFERENCES 


Armacost,  Robert  L,,  and  Fiacco,  Anthony  V,  1974.  Computational  experience 
in  sensitivity  analysis  for  nonlinear  programming.  Mathematical 
Programming . 6 t 301-326, 

, and  , 1975,  Second-order  parametric  sensitivity 

analysis  in  NLP  and  estimates  by  penalty  function  methods. 

Technical  Paper,  Serial  T-324,  The  Institute  for  Management 
Science  and  Engineering,  The  George  Washington  University. 


, and  , 1976.  NLP  sensitivity  for  R,H,S  perturbations  i 

A brief  survey  and  recent  second-order  extensions.  Technical  Paper, 
Serial  T-334.  The  Institute  for  Management  Science  and  Engineering, 
The  George  Washington  University, 

, and  Mylander,  W,  Charles,  1973.  A guide  to  a SUMT-Version  4 

computer  subroutine  for  implementing  sensitivity  analysis  in 
nonlinear  programming.  Technical  Paper,  Serial  T-287.  The  Institute 
for  Management  Science  and  Engineering,  The  George  Washington 
University, 

Colville,  A.  R,  1968.  A comparative  study  of  nonlinear  programming  codes. 
IBM  New  York  Scientific  Center  Technical  Report  320-2947. 

Fiacco,  Anthony  V.  1973.  Sensitivity  analysis  for  nonlinear  programming 
using  penalty  methods.  Technical  Paper,  Serial  T-275.  The 
Institute  for  Management  Science  and  Engineering,  The  George 
Washington  University, 

McCormick,  Garth  P.  1972.  Computational  aspects  of  nonlinear  programming 
solutions  to  large  scale  Inventory  problems.  Technical  Memorandum, 
Serial  TM-63488.  The  Institute  for  Management  Science  and 
Engineering,  The  George  Washington  University, 

Mylander,  W.  Charles,  Holmes,  Raymond  L, , and  McCormick,  Garth  P.  1971. 

A guide  to  SUMT- Version  4i  The  computer  program  implementing  the 
sequential  unconstrained  minimization  technique  for  nonlinear 
programming,  RAC- P-63,  Research  Analysis  Corporation, 

Schrady,  D,  A,,  and  Choe,  U.  C.  1971.  Models  for  multi-item  continuous 
review  inventory  policies  subject  to  constraints,  Naval  Research 
Logistics  Quarterly,  18i541-463. 

van  de  Panne,  C,,  and  Popp,  W,  1963.  Minimum-cost  cattle  feed  tinder 

probabilistic  protein  constraints.  Management  Science.  9i405-430. 


T-335 


COPY  AVAILABLE  TO  DDC  DOES  NOT 
PERMIT  FOLLY  LEGIBLE  PRODOCTION 


APPENDIX  A 

SUBROUTINES  SENS,  LMULT  AND  PRESEN 


CGCI 


035  2 
003 
0C'>* 
C33S 
OC  16 
05  5 * 

occo 


©CC9 
0015 
0*11 
C I » 
Of  I J 
<01* 
66  IS 
0019 
©5»y 

CC  J A 
C'J  V 
C325 
fdl 
oopp 

tfpj 

Off* 

C 3 2 A 

OtP» 
OOP* 
• OP* 

00  If 
•Oil 
ot  •/ 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

t 


c 


*U..n  JuT  l»*l  Sf  N$ 

15  19** 

This  W»:-SICN  *)*  TW  il  MMTIVITT  txm',11  SunftOUTtwF  IS  US*d  TO 
<fj»fWUTF  the  O J e ' c n t:**  L PI  A I VA  t I Vr  S (X  * tNO  r • | TM  *«rSP€CT  TO 
c • f * i u r %.  am.  ir*-  s cnn»  n in  t»»c  «»mt  paaipci.  thy  oi«‘CTIonal 

o*  » I V A T l V 1 % A.‘  'ST|MAT<0  # 0*  in*  KlktMU1-*  AT  A T l t*f  •IfM  NPAP  Of;  I NC 
Ti»r  »,u-  •.  t.  ,:f  Pic«wr  ims  ikAlvo  In  tm*  SENSITIVITY  analysis.  THE 
U'.»  '»*  Til*  PAPIM'-M  kS  PAP  Ip'  | MUST  »*  f.  HNS  I STf  NT  THT  OU'iMOUT  THC 

UV**»S  '-U'  - 'TUT  I N*  S. 

f»<r  su»n.uij*4»  is  unm>  *o«  a sinsitivity  analysis  at  thc  pinal  suo- 

P=r-LIP  O rrjp  * MnSITIvITy  analysis  at  tacm  SU6P*06L*n  alOnG  The 
M I N I p I / tN"«  T * A Jf  C TL">  Y.  Or*AC|2C»  IS  ▼ AAA  AY  0*  01 » * T M| NC I NG 
|lH<-Bv«lS  <-»«  r5P<jN^IN<,  »<»  T Mr.  PAR  AM*  f r r-S  PAM  20 1«  OPAW(ffl)  IS 
ASS  I C*  r U V *LUr  S |N  *.«>••*•  *:UMNr  P*«Ol»  . 

This  A.^r'AfH  TO  S*NSlT|V»Tf  ANALYSIS  IS  OUC  TO  Aa  V.  P)*CCO»  TH* 
p i *s  t v • s i jn  » *s  corro  hy  •*.  causct.  the  mcono  vipsion  «as  cooio 

OV  «.  C.  -YLAN1EA.  Tills  IS  ThC  T«|tO  VISION  «H|CH  |S  AN  laTtNSION 
r.¥  In’-  rr*.0  Vf  E »,  | r N t"  PHYIT  SENSITIVITY  ANALYSIS  ALONG  TMt 

MINIMIZIN'.  TPAj-CTrAY*  A*.P  «AS  COVO  MY  A.  L.  AAMACOST, 

IMPLICIT  a*  ala"! 

A*  AL*A  A NO  I V.BAt  |0.»  PS  I .TMr  TAf ,»rPl  . IfPf 

CLMVf  N/*HA<  /*»;»•  |,l.'  Lt?/M*AC?3.ni.N.«l.«<N.NP|  »NM| 

<t V/'OAL  ,**/ 

CCmmON/C«'TN  »/NT|  *M*  ,»»T  3.NTA.NT«.NTC,.NTT  tNT0,MTO*NT!0 

C *.*»'*  .*«/VAL  Ur  /f  ,E.,PE  | |VMO 

avavs/rfv?/(JLlM?l'  I .01  1 API  / 3 I.  A HOI  N.  PAT  tO.fPS  I * Th»  T SO* 

intc  iA  ,»  UM|*»|  ,iiizi*i.i?ipri,i!(;<>i,i«mf  i.soti**  i.p*i  » 

»»:  P.»*l  .r  I .o  J|  | AO  | ,l»OTT  A?«  J3  1 ,0IA6|  TO  It 

• V j , 4r»L«t»  S IGl  .Cl  ,NP».AS*  .ns  AT  is 
Ct,  .-*  S/S^N/P  AB  c 2CI  « I *P  4 A I t;  ) ,NPAB  , I ST  N$ 

CCY*ON/f  IPJPT/'.f  *0**l  . Nf  IOP?.  N*lOI*%NPll]P*|N(lOP|<tfT|»l(Pt 
IzImEnSMN  OUtUtliM  LTliniilpH|*OlillMiiei 
01  «»*.s:on  r'.LHUl  ACI.OUI  AO)  *«TCST<  251 
m »r.  »A  * l 
Nr  *h«s  * o 

IP!  S I - I .N 
OHUIII  • OIL*  I I I 
S 3'1’UI  • 0CL*M  I I 
<4|  L STH-f 
4**1  - l.C'-IO 

P**M#  I * • <| 

I V NS  a | ‘ *.NS  • I 
<AWL  PA.  Cl* 

A A I t * Out  I,  l<  ANO  PAP. 

..M'll.’CI  «H0 

|C  M'A»ATIIH  / / IS  » « 20**  ST  NSITIVITY  ANALYSIS  // 

1*.  A » 2>H  l»‘*  VALUE  C*  P|»**r,|  IS  atltaSI 

*r  rut  MAI  «/*.*.  A«M*<«r  POINT  AT  amiCm  1MC  KST|mat«  op  UNIIMVITT  AlLL 

I U(  -*U*  IS  I 

nu  a * !• i .n.a 

I 1«M|NSI I IS.AI 

VAIt'l...  101  IO.IUII,  ^Yl.ll) 

)T  PUAMATI  i(|l,)Y||,|)(|N|*,A|A,T|  | 

AC  tOST|M|* 

•A  lit  | , SO  I ill.  PAA  |t  l,|M*AA(|||.  I • I « AtPA*  ) 


016000 
oier.oc 
Cl AT^O 
016710 
C 16720 
Cl  67 30 
01*740 
CINTS3 
f I6YP0 
0167T0 
016700 
5 I 6*90 
0I60C0 
516810 
0ISO20 
01  AO  TO 
C I 6640 
016060 
C l(«8f,0 
01  66  7'J 
016690 
0 I 06*3 
0)6900 
016910 
016920 
016930 
010*40 
0I6AS0 
01*960 
016970 
010960 
0)6090 
0)7050 
OIPOIO 


01  PC  >6 
01*0  30 
I I *C  AO 
01 7050 
C I *.  66 
0 1 70 70 
01  Yf 63 
01 7090 
017166 
01*115 

omte 

01*1  3C 
0 I ft  AO 
01*160 
017100 
0171*0 
01*100 
<1*100 
•Dro 
01*0)0 
01 *520 


43 


COPY  AVAILABLE  TO  DOC  DOES  NOT 
PERMIT  FULLY  LEGIBLE  PRODUCTION 


T-335 


0031 

so  f ot  vat  t4*.M>  ptuMrtf*  vacuo  fl 

lIPPCfCNClNC  INTCKVAC  / 

0 1 T 2 SC 

1 « 1 * • 1 MS.4B.CI4 ,41  | 

• 

01*240 

o*n« 

4*.  1 T»  It  ,M» 

C 1*250 

co  is 

51  KI  **»  I//I 

01  * 2*0 

c evAcu»rr  r,  c ano  m. 

01*2*0 

C 0 16 

call  I'SINU  1.F  t 

* 

0I*29C 

CC.  1* 

CU  TO  *5 

01 *?«0 

(Cl' 

oo  a-  j*:  ,*»pn| 

eiTjro 

rc  j* 

CALL  A 1 SlMI^TjUII 

01*110 

00*.  *3 

•0  CCNTlsyr 

« 

Cl *320 

ce«t 

If  IV  BO*  ’•  * r 0 • 9 > CO  TO  OS 

0042 

rui  ppisemou.atcsti 

c CCWl’Ul*  0.  L #. 

01*930 

041 

05  LML  Ct-BCI  (0) 

01*940 

0044 

on  r,  i»i,n 

01*350 

0041 

aimiii  * oti  in 

01*900 

1046 

01  C'f.T  |4Ur 

01 *3*0 

C COIPUT*  (nr  ».•••*  P - STo»eo  In  A* 

01  * 3S0 

004* 

CAti.  I'-romi) 

01*393 

c PfKo-w  THf  c-u  OCCOWPOSIT ION  or  A. 

01*410 

0041 

OP  Ifl  IM.N 

0 I *4  » A 

CC4? 

0*  L ■(  1 l«C ,3 

01*420 

oov 

10*.  CO.f|NUr 

Cl*4)0 

fell 

MP.Nf 1 

f 1 *440 

C5'.2 

NT  I»1 

01*450 

CIO 

Calc  I NVl  AMI! 

C 1 *4  CO 

C C Hf  C 4 TO  «APt  SUP*  an  OPTHOCONAL  »OVC 

If  NOT  5TTCMPTCQ, 

01  *4  V 

0 0»-4 

on  lie  1 * i ,n 

01 *443 

r^ss 

If  (U1  L * 1 I I.CO.O.C)  CC  TO  113 

01*490 

OCCfc 

1 T'  < a .ICM 

f 1*510 

C35* 

1 45  ro»4l|t4i«,  THf  MBTbIK  Of  SCCONO  PAPTlALf  If  NOT 

POSIT  I VC  OCPINI.TC 

C»*S|0 

1.  T I VI  T V ANALYSIS  is  TEKMlNATCD  1 

C 175*3 

0*54 

CO  *•>  ?'  3? 

01*530 

0019 

lie  cr*.TiNt.f 

C 1 *5*0 

OP  l.s  1 »l  i NP A £ 

01*5S0 

COM 

B ?!•(  5 1 . MAC  i 1 1 

01*560 

r e»,2 

12C  CCo.lINi,* 

€1*5*3 

OCfcl 

0"  >f  1 J » I ,«*PAP 

01*550 

rci4 

|»  |M  Itc  .-*1  l.c  TO  ||S 

9 06*- 

|i  |*'r«,Tl  Jl,>  U.'  | CO  TO  201 

rose 

115  If  M.  « •A4.>  j.JI  CC  TO  121 

CCS* 

Call  t yul  T 1 1 .CILJMJ.O'  M.I  Ul 

< cona-ut r ot;»*L  PI/04IJ)  and  Dif|/0MJI 

UfINf  CCNTfAL 

OIPPCPlNClNft# 

01  *5SC 

Of  *4 

m PAaTjI  U lj|  • PPAK  J» 

< i*oco 

C CS» 

CALL  SMIMIJ.OI  » 

•1*5 10 

r cto 

lfl4.rc.CI  CC  TO  no 

Cl *t  2 0 

COM 

00  121  !'».*• 

01 *6  1C 

t?V 

call  4'  stnti | ,*J|  III 

01  *#46 

COM 

II  t A J( 1 t.CT. AVAL  1 CO  TO  120 

CI*6fC 

C4*4 

122  op**!  ji «'  .i  tr»AMji 

• 1 *960 

00*5 

PAHJI  » A2NIJI 

•1*5*0 

ec** 

• «-| If  <<  ,12*1  J.OPA4I 21 

01 *#40 

*S** 

124  fOMAUILH  WCS'TTINC  OPA»«  ,12. SMI* 

.*14, f» 

• 1*590 

CC  »• 

IMOPA^OI.CO.O.I  CO  TO  201 

017*00 

00»t 

CO  1 r 121 

•1**10 

0060 

125  COnTInuS 

•1**20 

0'J#1 

120  I2IW2.IO.0I  CO  TO  t»ft 

•1**10 

-00# 

OO  |2T  1*1* *1 

•1**40 

et4» 

W/P|*W4 1 

•1**93 

COM 

CALL  •*  51NT|W#P|  «»J|W|P|  1 1 

•1 **50 

0 3 *4 

12 Y CONflNUl 

•»***• 

cc  1* 

It*  |»  (N»  411x4,1  O.f.l  CO  TO  120 

0C4* 

CALL  L*U  *11  .0»L*U.P*w. tHf» 

OM 

It*  C*tL  C-ACI2I 

(Ml 

00  1*0  |4».N 

• 1***0 

ee>: 

Oil  #1  1 )•'>!  L *SI  1 1 

•1*000 

C09I 

113  CUNTINU* 

01 *010 

0012 

Of  mMtPPiflJI 

•1*090 

0441 

PAfc  ( iliru  1 J | - 6(4 

C * »•  *4 

06  » 

CALL  *rSTNMS.VALl 

•1 *040 

ten 

|f|w,»g.f»  CO  1Q  Ilf 

C 1 *050 

COM 

on  1*5  I*|.w 

0»  *05T 

CCS* 

(ALL  a*  SlNM  t.AJI  III 

•1 *0*0 

00^4 

If (*2t 1 I.CT.AVALl  CO  to  Ilf 

•1 *0«3 

?e»» 

CO  TO  12* 

• 1*966 

• 1 10 

US  CCN?|NOf 

• 1 *9'9 

0191 

|H  Mlw2.fU.SI  CO  TO  IIS 

•1*010 

0192 

DO  19*  1*1 .«! 

• 1 *020 

0101 
01  44 

4|2|*4. | 

• 1*0 1C 
•I*f4« 

OIOS 

|1V  (P4l|4.»* 

•1*05* 

- UU  - 


W AVAILABLE  TO  DDC  DOES  NOT,. 
PEBMIT  TOUT  LEGIBLE  PBOOOCTION 


T-335 


0106 

ISO 

IFI4'.  «CI  4.*0'  *1  CC  ?0  ISO 

C107 

C*u  l **VL.  til* l>t  L •'U  » 04  M . pu  | 

Cl  06 

1 SO 

CHL  c«  »ri?» 

01  'i*l 

i.»« I r>r  <w*i  i/t)rm 

017*70 

Cl  1C 

M i»:  m.n 

017040 

01  1 1 

u'hiii'Ii.'lhii  - ot  l*oi  1 1 i/ocm 

©170)0 

0111 

140 

Ci  At  I4ur 

ODC^O 

01  1 1 

*>*•  1 Jl  >•  «>I<U» 

OIOC  10 

c w»im,  *in  *nr  or  iiimo  4*  r.»n.  vf  a • i>#  ro<*  *» 

01  AC  10 

C »M  t.r  h'U'L*  *40  Jl, 

oi  oeso 

4114 

C 41,  c I'.vl'MAI 

C 1 40*0 

C *4 INT  nut  t»/f  *1)1 

01 0040 

0114 

• HtM'.,IS(l  J 

01*060 

0114 

140 

»fj>-M«M«7M  i-firMtAiivri  *•<  • 1 1 m »c§*ccf 

OIOC  70 

01  1 7 

00  17*  l«t.N.4 

010000 

0114 

1 1 t«lNM  t » *.N) 

OIOC 00 

"1  14 

• 1 J J.O’ l«f  JJI  . JJ*l.ttl 

ototro 

CI14 

106 

P'*f7»t<  # 1 4 m Oil  . IP  «/•<>  *iM  4«7|  1 

01*110 

cm 

l»C 

CCST|<U». 

010110 

0111 

iriNf irr«,fo,ei  CO  to  m 

ot  ?y 

IM*.*  (,.»  » '.c  t r ir*! 

Cl  14 

Cut.  •.•'ui,  1 1 t.oi’t  *u.ct  **.r»u» 

C 1 14 

■ » 1 t*  1*.  it-?*  j 

C 114 

)40 

POM»M/  * 1 **  U-0fB|V4f|Vf  ft  0 1 t M HtTfCt 

TO  M»i»CT(l  • If! 

0117 

00  1* 1 1*1 

9 

0114 

II  * n|n'I 1 *4.*» 

C110 

• ■»»'«»  . *'?)  ll))*t>u(J)II.JJ*l»l|| 

01  15 

S4» 

Out  • 1 1.!*!*.  Cl  *«7  |» 

01  11 

SSI 

cr».»  i nu* 

01  11 

SOI 

l«  <«<*rC.C 1 CO  IC  17*  . 

nil 

C4|  t CULM*  •0<L*V*C*M,0U> 

01) 

itnMt.mi  j 

, ^ 

oil'. 

54  0 

4 | *•  ••OI*l«lt|VtS  «|tN  i||A«t 

TO  MiAMfllt  lilt 

Cl’* 

bU  )t  1 1 * 1 • •*/  .6 

e i ■ 7 

it  • !•»•-*<> 

Cl  >4 

• Ml'  »•  . AM  i t JJ.'JUi  JJ**n  • J)«t.tll 

ei  i* 

SOI 

Pit  »n(<Uf  0*1  . IZ.lM)  4.014. 711 

ci*: 

sot 

t l'«o* 

0141 

S 74 

cn*.»  jsur 

C o»/o*ui. 

01 01 13 

C 1 41 

or  in  i*i. v 

ct*i*e 

01*1 

or  • jr  • «!»•(  II 407X11 11 

0|*4 

10' 

<C'*t|Mj» 

ei*i*3 

0 7*  I«|t  Uf/t  *l5l. 

Cl  M *0 

C l *4 

* 

O' 

* 

01*1 *3 

31*4 

IV* 

rci**M/  pi,i)Hnrinti)/r*»  .ci*«0/10m  •< 

©|4»v 

C 1 4 7 

IOC 

C t *•  1 1 * 

«i*»r.f 

Cl*4 

60  to  7*1 

•14710 

014) 

101 

•*itrt*.ic*i  j 

nut; 

4l4t 

104 

014110 

3141 

40  tO  IOC 

0410 

0141 

101 

NT  •••*»» 

CI014C 

014) 

C*IL  ••  )7  C 7 

OIOK  3 

Cl  44 

OO  10*.  tal.N 

<1*1*0 

eif. 

OUtlll  * I 1 

<101*3 

0144 

104 

OtC*  »l  1 1 • OUMI  1 

Cl  410 

0147 

4*7  TU*N 

•too 

0144 

• NO 

01031) 

T-335 


pnpY  IVIIUBIE  TO  DOG  DOES  NO] 

Sit foil?  legik  p«cW 


CO-)! 

lVl*  JUl  |N*  t*ULT|  IKO.O*  L«U.Q«V.OVI 

c o i ■» 

l«Ml(|!  »**L«*f  *-»**'W| 

COC  i 

ci  *c»*  khCh.sMic.rm  .fmo*c  .if pi  « «r*»* 

OCi« 

CU*‘,*v',i/V*»‘l/u^  * .r-M  i ?->».»< . »“*»  lfcl«,miiN*liN«| 

oe*r. 

(r*»«i  I't/KHuf  ft  .o.p*'  .«  lir,M«,ijKO  1 #n«0 

OOOC 

tOmm  .N/tOtl/  h,  mi.  «•/ 

oev 

CC*»*.v  l/C6  0f/0»CM9fl  ,U»  t*'  t?C  1.0*0 IN. .?HCT*0 

.M  .f  1 .»  JM*->I  ,»*rTT,*H.6*(M99l  .01  *C< 99). 

JPi  * v'.4l)'L*>41IC|.M.  .NO*  T 10 

?o*>« 

9 1 t*rh't  t On 

C9  JO 

m*i  i w » w t i 

Of  n 

60  TO  ( t .9. 1.1) . INO 

C INO 

• r 

on 

00  199  |a|.»"9 

001  * 

l^O 

f'l  l «w<  II  • 9. 

OMJ 

Cf  1U*N 

C fND 

■ t 

691* 

1 

or*  O',  in.***/ 

COI*!. 

oo 

0'  l *>0(  1 » » H J(  | 1 

0616 

M.  TO^N 

C INO 

• 9 

con 

9 

1)0  1 •)  .***/ 

;tu 

ec 

Ud.“jtl)  • (O'LMUtll  - •Jtlll/OCH 

CCIM 

M.1W-  h 

0C99 

s 

DO  yr  lat.WM 1 

901 

(*u  6u*ntti ) 

c 09 

C*CC  f t5»NT| I.VCll 

CO  9) 

SUM  * 0. 

:ei« 

C9  »l  99*).N 

o:»? 

*1 

*o«*  » 00**  • DCUilianciion 

(CH 

if  i iNn.ro. *i  6C  ic  •*' 

rc*t 

owl  n • 6 1 "W»  1 » -**i*aft) 

Cf  l» 

66  fc  T9 

JOM 

• 9 

OUtll  «*  ?.*ISU>  • Of  L«Ut  1 1 1 /•»*♦) 

00  >0 

TO 

CCNTlNUf 

00  ft 

• 7 »g*>t 

00  19 

ifcO 

46 


non 

<U''(OU*  IN'  Pf>r  "Vi  Mfiu.xTFSl  1 

cc  •> 2 

IMPLICIT  Ci al'MC A-H.L-Cl 

C C 1 3 

C '.mm.-j/'  .uh./<|,’i  1 .filLC  23  1 » A C 2'  » 2C1  .N.M.MN.NP1  .NHt 

000* 

cr**  .h/w*lub /r  rr, .t  s igma.r  jc*o » .»mo 

coo* 

(J/M,/  It,  H|  ,1*1 

0036 

CCMMCN/M  n/p*h  C 20  1 , DM  *t-  ( 20  1 .NPA*.  ISC  NS 

or  t r 

D!  >*'  ‘.5  l r>f.  r, * ( *c  I .OUl*OI  .K TEST! 20  1 .KLISMZO) 

et  36 

MPM  2 • M • M/ 

0"33 

FT' 61  « C.'Jl  • SAMS  IF  1 

C 0 1 3 

00  1*9  J*|.NP*ft 

0911 

AUSTCJ1  « J 

OC  12 

PAMJI  p*MJ>  ♦ OP  Aft  ( J ) 

0013 

CAl  L ei.STNTIO  ,07  1 

OCI* 

|r(MPM/,rc.?l  CO  TO  20 

oms 

00  I*  1*1. MPM* 

0616 

C»IL  cr  SI N T ( 1 ,OUC  1 1 1 

oo'  t 

1C 

cohMNur  # 

06  1 8 

20 

0F^»  * 2.  • OP  Aft  ( j | 

r.i« 

P*M  Jl  r Pta(j)  - OF» 

0026 

CAl L ftrSTNTCO.MFI 

9021 

ir 0.31  CC  TO  AO 

C 322 

00  33  1 • 1 .MP-J 

CC23 

cai  i_  c'-srNM  i .exit  1 1 

0C2* 

30 

com  i Ny 

0026 

*9 

OFfps  ■ 1 OF  - xfi/ofm 

0026 

IF  I mom/,'  0.3  1 CO  TO  60 

0C*7 

OO  SO  1 «l  .MP*»* 

0020 

50 

OUII1  ' C0UC11  - GKCIH/OFM 

0029 

60 

SU*  « CFtps 

OC  33 

IFlv.t  0.01  CC  TO  80 

00  31 

OC  t'i  t»l.V 

CC  32 

7C 

sum  * sum  _ BMC/p jc i i*ouc n 

CO  33 

80 

IFCM2.ro. 01  CO  TO  95 

OC  3* 

TSUM  * 0. 

or  36 

on  l«i .“2 

C9  36 

1 M « I «M 

00  IT 

90 

TSUM  * TSUM  ♦ r JC  |M|*OUC IM) 

OC  33 

SUM  SUM  ♦ TSUM  • ? • /ft HO 

0 0 39 

96 

OfLCJI  * SUM 

06*9 

PA-JCJI  * PAPCJ1  ♦ OP  A®  C J 1 

09*1 

DT'ST  * CAOSCOrtCJI) 

00*2 

If  COTFSl.C'.FTFSTI  KTCSTUI  m | 

OC  * 3 

109 

CUM  INUL 

00** 

• ft  H*  C«  «<.  V I 

99*5 

600 

rO'-MATC  22*  . 3*M9PT1MAL  VALUE  FUNCTION  SENSITIVITY  // 1 

9C*6 

D'J  259  | * j ,np*.  ,5 

33*7 

| | *M|*1'  » | «4  ,SD*C  1 

C0*8 

»R|l?C6.'CI  1 C C JJ.OELC JJ1 1 . JjMt.ltl 

0C*9 

601 

ro»-M*Tl  SCTH  or /OAC  • t2.2»«l«.61**7)  | 

OC  V. 

too 

CONT  JNU' 

OCSI 

JJ  • 0 

0062 

OO  2S0  J«1 . NPAP 

CCf  3 

IFCKTC  STC  Jl  .F.O.C)  CO  TO  250 

©C5* 

JJ  • JJ  • 1 

C066 

RII5HJJ)  * J 

0068 

250 

continue 

9667 

1FCJJ.C6. 0)  CO  TO  300 

0066 

• 6 1 1C  C 6 • 602  I * 

1 060 

*02 

r0CM*T(  /SI  h CFTAILIO  SCNSIHVITV  FCSULTS  FQU.O*  FO*  FA*ANCTC*S 

toe: 

• ft|ircs.603l  1 *l  isle l|.  I-I.JJI 

06  61 

60  3 

FU'-maICIh  4‘jI  I3.2M  ,)) 

0062 

• »l  l*  ll.'Ctl 

0063 

*0* 

FOAMAIC/I 

006* 

M T U*  N 

6066 

300 

»*  I vr  c 6 • 1 05 1 

9966 

66  5 

FOftMArc  « .*•<  TMtVF  l«F  NO  OFTAILCO  *«N6IT|V|TV  MULIt 

//» 

606? 

•ElUAN 

0666 

f NO 

1*7 


T-335 


i 


APPENDIX  B 


USER  SUBROUTINES  FOR  SCHRADT-CHOE  PROBLQt 


com 

SliP*OUT|Nf  KCAD1N 

CC32 

ILLICIT  pr AC**< A-h.Q-l | 

0C03 

cn^M'ix/isy/  ir  ta(?oi  iPhk?'ii  «nrNSE(20  i«  iocntoakhi 

GC34 

Ct>M*0N/SSN/»»A3|  ?C  1 ,DP*P  <201  ,NP„B  , 1S€N* 

005 

901 

FCPM*T(15M^12.C  » 

0006 

BC*0<S.9CII  M .PACU  1 ,PAP(2I 

OOOT 

«(l  Tr  | Nl.PtPI  1 ),PAQ(2I 

0006 

0(1  l-*>5  1 «l  .Nl 

0009 

Sf /0( 5,901 1 lOf  NTl‘|  ).  (PAP  (4*1-2*31  , 3»l  ,*> 

0010 

i,Pt  tci  5,901  1 (Ot  NT  Ml,  (r»#j«»|-?tj|,j«l,*| 

0011 

too 

CCNT INUF 

0012 

NPAF  ■ t«NI«2 

0013 

PCTUFN 

COM 

FNO 

coot 

SuePUOTlNF  PFSTNfMN.VAtl 

00  02 

INOUCIT  ft  ACM ( »-n,n-21 

C?03 

COO* 

(C‘«if)N/5H»fir/x(?r  1 , DF  L ( 23  1 , * ( 20  , ?C  1 , N.O.  NN.NPl  ,NM| 

'‘005 

CCNM'  N/tNV/u'TAl  >0  1 .PHI  <20  I ,DLN$E<29  I.  IOCNT(t0l*NI 

000* 

CC^ON/S'N/PAM  2C  | ,OP*P(20I,NP*P,  ISFNS 

COOT 

v*c  ■ o* 

C95* 

IF  1 IN. f 0,01  GO  TO  300 

oro<* 

IF  1 IN, f 0.1 1 CC  TO  100 

C010 

2C0 

00  250  1*1. Nl 

cou 

13  ■ 2*1-1 

0012 

OO  ■ K<  1 31 

0013 

IF  (00,1.1  .9.1  GO  TO  1 80 

CCl* 

230 

VAC  * V*C  ♦ PAB <**1*21/00 

0C15 

V*L  • P*»l?»  - VAC 

cot* 

etTUBN 

00  1 T 

MO 

VAC  ■ -1 .0 

COM 

*t  3U*N 

0019 

100 

00  1 55  l*l.N| 

0020 

131  ■ 2*1 

0021 

1 J m |Jt  - | 

0022 

«(  * MlJtl 

Of  21 

IF  (o«>  .1  r ,o.  1 CO  TO  MO 

002* 

00  a Kill 

CC25 

irtOO.UT.O.  1 GO  TC  MO 

002* 

160 

V<L  • VH  ♦ p*p<**I*I  1* < •* *00/2, - P** <*•!•!  1 1 

0022 

V*l  • PtKII  - VAC 

002* 

Of  3U*N 

f 0 29 

390 

00  J'O  l*|.  Nl 

0939 

131  • 2*1 

0031 

13  • 131  - 1 

CO  32 

00  • >1131 

CO  33 

•■•*<1311 

00  3* 
0035 

ug  • PAP (*• |«1 1 

003* 
CC  IT 

53  a P AC <« a | | 

CC  *5 
fl(  13 

OF  C T*  • *P  - uU 

00*0 

IN  • OF  1 TA  / 33 

00*1 

C*LC  ASOTfM/N.r  1 .UtM 

CC*2 

PM||  1 1 •.  F | 

00*3 

UCN  .F  m • 01  N 

09** 

Of  TA(  1 | a*.  5*1  M3*55*OrtTA*oec2*l*#H||  | l-33*0«CT  KHNll  Ml  1 

09*9 

330 

V*C  a VAC  • or  TA< | 1/00 

fO*t 

Of  than 

00*2 

CNO 

1 

] 


i 


- 48  - 


T-335 


00)1 

c::2 

ccoi 

OC)* 

oco6 

co)6 

0332 

coon 
C 039 
0C1  3 
001  I 
0612 
CM? 
0014 
CO  I 5 
COl* 
001  2 
0 014 
OCI) 
0020 
0021 
0022 
0C2) 
0024 
CO  20 
C 0 26 
0C2T 
C4?4 
0C20 
00)3 
cost 

CO  12 
CO) 
001* 
C 0 35 
0326 
OC  3» 
C 0)0 
00)0 
0040 


0CC1 
CC  02 
CC  33 
0034 
0003 
00  06 
00)2 
00)4 
/'(')O 
©Cl) 
0011 
0012 

0013 

0014 
C 0 I 6 
OC  16 
001  2 
coin 

OC  !*• 
0026 
ceil 
C)?2 
C C 23 
0C24 
CO  25 
OC  16 
022 
CC?) 
CC20 
OC  J) 
00)1 
00  3? 
05)3 
00)4 
0033 


CO)l 

oeot 

003) 

0004 
0003 
0 006 

OOOT 

0006 

eooo 

eoio 


sun^ounsf  06401  uni 

IMPLICIT  OC*L*X|A-M.O-Z> 

CCkm  ;N/<5M»cr.  /x|*f  » .02  LI  20)  . 4 1 20.20)  .N.M.MN.NPI  ,NM| 
rn**M,  N/tsv/  ♦ *.  f 4 1 ?o ) . ph 1 1 20 ) .oeN$ei?o  )•  iocntiio  >.ni 

CC*"  lU/'.f  N/P *W<  ) .DP 4P| 20) ,NPAS  * I SENS 

IM lN.fO.0)  00  TC  300 
|F|  IN.FU.l  ) CO  TO  100 
DO  236  1*1 .N! 

I J I * ? 4 | 

13  * 1 J I - 1 
00  * XI t J> 

OC L I l Jl  ) • 0. 

250  0‘LIIJ)  • P*w<4*  (*)  1/00/00 

R6TU6  ’4 

100  DO  150  1*1. N| 

I J 1 * 2*1 
IJ  « Ul  - I 
on  i iji  ) t -p«b  («•].{  j 
150  orti  I J)  » otn  I J ) 1/2. 
f *•  Tyus 

300  00  163  I « I , N I 

121  * 2*1 
IJ  « IJI  - 1 
00  ■ XI  I J) 

PP  * XI IJI  1 

UU  • PAP | 4*1-1 » 

35  ■ PAP|4«|  ) 

0Flta  * an  - uu 
2N  « Ot  L T A / 55 
CALL  ANOTPI/N.PI .OfN) 

PHI  I I ) a F| 

Ot  NV  I I ) ■ OCN 

OFTai 1 ) *0.5«u  S5«S5*Of LTA*OElTA|*PM|( I )-S540€LTA*DCNSE< II) 
Oft  <1 Jl )»<D1LTA*PM1< t |«  55  *0* NSF I f ) 1/00 

33ff  on. 1 1 j ) » -orTAi  i ) /oo/oo 

PF.TU«-N 

END 


SunonuTjvr  **Tb  | a(  | n.  |xx  » 

implicit  p» *l«p( a-m.o-i i 

CO**  1N/$MAVF/X(?'<t  tfLI  20I.AI?3.2)).Nf«I.HN.NP|  .NM) 

COMMON/ | NV/ 11*2  4 I ?r>|.  PH  1120  1.  OtNSf.t  23  I.  IOENT120I.NI 

CO*'M'JN/SrN/PAA  120)  .DPAPI2C)  .NPAR.ISCNO 

1 A I t .*•  « L O « 0 ) CO  TO  300 

tFUN.10.1)  CO  TC  100 

00  250  I4I.NI 

IJ  - ?*l-| 

00  ■ XI  I J) 

250  AIIJ.IJ)  • -2.4P4P(4«! 421/004*3 
A4TUTN 
1)0  I x x ■ 1 
6 C T UP  N 

300  CO  4 ,0  1 * t . N I 
IJI  » 2 * I 
13  • IJI  - l 
OO  ■ X I | Jl 
4-A  ■ X|  | J|  ) 

UU  * PA* I 4 • | • | ) 

33  ■ P AP 14*1) 

Of  C T A ■ I.C  - UU 
2N  • 0»  LT A / 56 
c*il  ammm  ja.f  i ,rrN» 
mh ) i i i • r i 
OF.N5*  II)  ■ OCA 

OFTAI  I ) •4.6*1  I S3*S6»OfLTH40fLTA) *PH| I 1 I • 33*0f L T A40CNSC III) 
AIIJ.IJ)  • ?..!♦»  »A|  | 1/00443 

AIIJ.IJ)  )•- 1(51.1  »».PM|  ( | ).S3*oeN3t(t  ) 1/00/00 
330  All j) .Iji ) • 0*1(11/03 
Pf  tyf.N 
tNO 


3UM0^gt|...  tsu  »•  I x x , Ph|  «U4  N3K  | 

IMPLICIT  «'»  *L»P(  a-h.c-21 
AX  • 0AX6IXX) 

T • I .3/1 1 • -»»•'*  2 3164*1  0*  A * I 
OI.NSI  • I . *P< -X««XX/2. 1 

PHI  . (if  nm  • 1 • I I 1 1 1 • 3 102 24 • f >t .621 256 )4T«1  ,r«l 4 261  *T 
1 - e.JM-.M4l*T  • 0.3)6)6131 
IFIXXI  1.2.2 

1 PH|  • 1.0  - PHI 

2 PftUFN  # 

«NO 


49  - 


T-335 


THE  GEORGE  WASHINGTON  UNIVERSITY 
Program  in  Logistics 
Distribution  List  for  Technical  Papers 


The  George  Washington  University 
Office  of  Sponsored  Research 
Library 

Vice  President  H.  F.  Bright 
l)ean  Harold  Liebowitz 
Mr  J.  Frank  Douhlcday 

ONR 

Chief  of  Naval  Research  (Code  430D) 
Contract  Administrator,  SLA 

OPNAV 

OP-40 

DCNO,  Logistics 
Navy  Dept  Library 
OP-9 1 1 
OP-964 

Naval  Applied  Science  Lab 

Naval  Aviation  Integrated  Log  Support 

NAVCOSSACT 

Naval  Cmd  SysSup  Activity  Tech  Library 

Naval  Electronics  Lab  Library 

Naval  Facilities  Eng  Cmd  Tech  Library 

Naval  Ordnance  Station 
Louisville,  Ky. 

Indian  Head.  Md. 

Naval  Ordnance  Sys  Cmd  Library 

Naval  Research  Branch  Office 
Boston 
Chicago 
New  York 
Pasadena 
San  Francisco 

Naval  Research  Lab 
Tech  Info  Div 

Library.  Code  2029  (QNKL) 

Naval  Ship  Engng  Center 
Philadelphia.  Pa. 

Hyattsville.  Md. 

Naval  Ship  Res  & Dev  Center 

Naval  Sea  Systems  Command 
Tech  Library 
Code  073 

Naval  Supply  Systems  Command 
Library 

Capt  W.  T.  Nash 

Naval  War  College  Library 
Newport 

SUPERS  Tech  Library 
FMSO 

Integrated  Sea  Lift  Study 
USN  Ammo  Depot  Earle 

USN  Postgrad  School  Monterey 
Library 

Dr  Jack  R.  Bursting 
Prof  C.  R.  Jones 

HQ.  USMC 
Commandant.  USMC 

Marine  Corps  .^houl  Quant ico 
Landing  Force  Dev  Ctr 
Logistics  Officer 

Armed  Forces  Industrial  College 
Armed  Forces  Staff  College 

Army  War  College  Library 
Carlisle  Barracks 

Army  Cmd  & (ien  Staff  College 

HQ.  US  Army  Trans  Mat  CommanJ 
TCMAC,  ADST 


Army  Logistics  Mgt  Center 
Fort  Lee 

US  Army  Electronic  Command 
Philadelphia 

Commanding  Officer,  USALDSRA 

New  Cumberland  Army  Depot 

F'rankford  Arsenal 

HQ.  USAF 

AFADS-3 

Maxwell  APB  Library 

Wright-Patterson  A I B 

HQ,  AF  Log  Command 
Research  Sch  Log 

Defense  Documentation  Center 

National  Bureau  of  Standards 
Dr  FI.  W.  Cannon 
Dr  Joan  Rosenblatt 

National  Science  Foundation 

National  Security  Agency 

WSF.G 

British  Navy  Staff 

Logistics.  OK  Analysis  Establishment 

National  Defense  Hdqtrs.  Ottawa 

American  Power  Jet  Co 

George  Chernowitz 

ARCON  Corp 

General  Dynamics.  Pomona 

General  Research  Corp 
Dr  Hugh  Cole 
Library 

PUnnuTg  Research  Corp 
Los  Angeles 

Rand  Corporation 

Logistics  Dept. 

Library 

Carnegie- Mellon  University 
Dean  H.  A.  Simon 
Prof  G.  Thompson 

Case  Western  Reserve  University 
Prof  B.  V.  Dean 
Prof  John  R.  Isbell 
Prof  M Mesarovic 
ProfS.  Zacks 

Cornell  University 

Prof  R.  E.  Bechhofer 
Prof  R.  W.  Conway 
Prof  J.  Kiefer 
Prof  Andrew  Schultz.  Jr. 

Cowles  Foundation  for  Research 
Library 

Prof  Herbert  Scarf 
Prof  Martin  Shubik 

Florida  State  University 
Prof  R.  A.  Bradley 
Prof  I.  R.  Savage 

Harvard  University 

Prof  K.  J.  Arrow 
Prof  W.  G.  Cochran 
Prof  Arthur  Schleifer.  Jr. 

New  York  University 

Prof  O.  Morgenslern 

Princeton  Un^ersity 

Pro'  W.  Tucket 
Pro  I.  W.  Tukey 
Prol  Geoffrey  S Wjtion 


Purdue  University 

Prof  S.  S.  Gupt i 

Prof  H.  Rubin 

Prof  Andrew  Whinston 

Stanford  University 

Prof  T.  W.  Anderson 
Prof  H.  Chernoff 
Prof  G.  B.  Dantzi g 
Prof  F.  S.  Htlher 
Prof  D.  L.  Iglehart 
Prof  Samuel  Karlin 
Prof  G.  J.  Liebcrman 
Prof  Herbert  Solomon 
Prof  A.  F.  Veinott,  Jr. 

University  of  California,  Berkeley 
Prof  R.  E.  Barlow 
Prof  0.  Gale 

Prof  Rosedith  Sitgreaves 
Prof  L.  M.  Tichvinsky 

University  of  California,  Los  Angeles 
Prof  J.  R.  Jackson 
Prof  Jacob  Marschak 
Prof  R.  R.  O Neill 
Numerical  Analysis  Res  Librarian 

University  of  North  Carolina 
Prof  W.  L.  Smith 
Prof  M.  R.  Leadbetter 

University  of  Pennsylvania 
Prof  Russell  Ackoff 
Prof  Thomas  L.  Saaty 

University  of  Texas 

Prof  A.  Charnes 
Prof  A.  C.  Stedry 

Yale  University 

Prof  F.  J.  An  scorn  be 
Dept  of  Admm  Sciences 

Prof  Z.  W.  Birnbaum 
University  of  Washington 

Prof  B.  H.  Bissinger 

The  Pennsylvania  State  University 

Prof  Seth  Bonder 
University  of  Michigan 

Prof  G.  E.  P.  Box 
University  of  Wisconsin 

Dr  Jerome  Bracken 
Institute  for  Defense  Analyses 

Mr  Wallace  M.  Cohen 
US  General  Accounting  Office 

Prof  C.  Herman 
Columbia  University 

Prof  H.  F.  Dodge 
Rutgers-The  State  University 

Prof  Paul  S.  Dwyer 
Bradenton,  H. 

Prof  J.  Gani 

University  of  Sheffield,  England 

Or  Saul  I.  Gass 
Mathematica 

Dr  Donald  P.  Gaver 
Carmel,  California 


Prof  J.  F.  Hannan 
Michigan  State  University 

Prof  H.  O.  Hartley 
Texas  A & M Foundation 

Mr  Gerald  K.  Hein 

NASA,  Lewis  Research  Center 

Prof  W.  M.  Hirsch 
Courant  Institute 

Dr  Alan  J.  Hoffman 
IBM,  Yorktown  Heights 

Dr  Rudolf  Husser 
University  of  Bern,  Switzerland 

Prof  J.  H.  K.  Kao 

Polytech  Institute  of  New  York 

Prof  W.  Kruskal 
University  of  Chicago 

ProfC.  E.  Lemke 
Rensselaer  Polytech  Institute 

Prof  Steven  Nahmus 
University  of  Pittsburgh 

Library 

National  Academy  of  Sciences 

Prof  D.  B.  Owen 

Southern  Methodist  University 

Prof  E.  Parzen 

State  University  New  York,  Buffalo 

Prof  H.  O.  Posten 
University  of  Connecticut 

Prof  R.  Remage,  Jr. 

University  of  Delaware 

Dr  Fred  Rigby 
Texas  Tech  College 

Mr  David  Rosenblatt 
Washington,  D.  C. 

Prof  M.  Rosenblatt 

University  of  California,  San  Diego 

Prof  Alan  J.  Rowe 

University  of  Southern  California 

Prof  A.  H.  Rubenstein 
Northwestern  University 

Dr  M.  E.  Salveson 
West  Los  Angeles 

Mr  N.  Sayal 
Silver  Spring 

Prof  R.  M.  Thrall 
Rice  University 

Dr  S.  Vajda 

University  of  Birmingham.  England 

Prof  T.  M.  Whitin 
Wesleyan  University 

Mr  Marshall  K Wood 
National  Planning  Association 

Prof  Max  A.  Woodbury 
Duke  University 

Prof  Jacob  Wolfowitz 
University  of  Illinois 


December  1974 


THE  GEORGE  WASHINGTON  UNIVERSITY 


,r 


BENEATH  THIS  PLAQUE 
^ -v.  IS  BURIED 

a^ Vault  for  the  iuture 

IN  THE  YEAR  20*36 


%«•«*  tl 


*V  4 


I#, 


v» 


TJTE  STORY  OF  ENGINEfHING  IN  I H I J>  YLA'K  OF  IMF.  PLACING  OF  IHh  V AVI  Li  AND 
ENGINEERING  HOPLS  FOR  THE  TOMORROWS  AS  WKH'ILN  IN  THE  REGORUS  OF  THE 
FOLLOWING  GOVERNMENTAL  AND  P HOF  FISS  ION  A L ENGINEERING  ORGAN  1 2:  A I IONS  ANt* 
/THOSE  OF  THIS  GEORGE  WASHINGTON  UNIVERSITY. 


r* 


, . BOARD  OF  COMMISSIONERS  DISTRICT  OF  COLUMBIA 

-UNITED  STATES  ATOMIC  ENE  R CjY  COMMISSION  ✓ 

DEPARTMENT  OF  THE  ARMY  UNTTED  STATES  OF  AMERICA  . / I 

DEPARTMENT  OF  THE  NAVY.  UNITED  STATES  OF  AMLRlCA  . ^ 

DEPARTMENT  OF  THE  APR  FORCE  UNITED  STATES  OF  AMERICA 
NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUT  ICS  : 

NATIONAL  BUREAU  OF  STANDARDS  U S DEPARTMENT  OF  COMMERCES 
AMERICAN  SOCIETY  OF  CIVIL  ENGINE  E"R  S 
, AMERICAN  INSTITUTE  OF  ELECTRICAL  ENGINEERS 
i THE  AMERICAN  SOCIETY  OF  MECHANICAL  ENGINEERS 
* THE  SOCIETY  OF  AMERICAN  MILITARY  ENGINEERS 

AMERICAN  INSTITUTE  OF  MINING  & ’METALLURGICAL  ENGINEERS  . 
DISTRICT  OE  COLUMBIA  SOCIETY  OF  PROFESSIONAL  E N G l N EERS^I  Nc  ’* 

THE  INSTITUTE  OF  RADIO  ENGINEERS  INC  ^ 

THE  CHEMICAL  ENGINEERS  CLUB  OF  WASHINGTON  , 

WASHINCTON  SOCIETY  OF  ENGINEERS  tWmmm 

FAClKNER  KJNC  SBUkY  a STENHOUSE  * ARCH1  TECTS^ 

CHARLES  H TOMPKINS  COM  PAN  Y - Mil  Lbtf&S  < ~ '%■ 

socttay  Of  WOMEN  engineers  * ^ -t 

NATIONAL  ACADEMY  OF  SCIENCES,  NATIONAL  " 

■'  / - 1S2*  4 r " >«cv*4^ 


Vi.  I 


TV 


a.*  T r 

m,' 


RCH  .cm:  Sens*  ' W., 

THE  purpose  of  this  vault  IS  lfgp:i[mtT.,&^  ^t>  PS  DEDTCATBD.  ft* 
CHARLES  HOOK  , T oTT. R K^l  N S?. 6 oWirfr.  or  ft N-cfflEtm t\c?. 


JCGAUSe 
tDTlM  I j NT 

*r~ 

iS^scntKpri 

m * -.si 

T* 


glflS  ENGINEERING  COrmfiTOMl 
HIS  NATION  AND  ToS®flt|||S 

BY  THE  GEORgV.V^HiN'i 

nc  • . 

board  of  r r t j s t pt * Yv 

I U S F TH'  » 4w 

v - nv^ 


fc  • ■»;*  '*»■■■  •*?  *-•'  3J?  «P 

■ ■ " -■■  fc  . > :*->  % s* 


\ -l  T 


fA:  V V- vSW  •*:*.-* 


•;,V 

->*sdMas 


• *t  I 


' I 


To  cope  with  the  expanding  technology,  our  society  must 
be  assured  of  a continuing  supply  of  rigorously  trained 
and  educ  ated  engineers.  I he  School  ol  Engineering  and 
Applied  Science  is  completely  committed  to  this  ob- 
jective. 


